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1 Chapter 1 - Introduction 
1.1 Motivation 

This report, concluding a one-year NIAC Phase I study, describes a new structural and 
mechanical technique aimed at reducing the mass and increasing the deployed-to-stowed 
length and volume ratios of spacecraft systems. This technique uses the magnetic fields 
generated by electrical current passing through coils of high-temperature superconductors 
(HTSs) to support spacecraft structures and deploy them to operational configurations from 
their stowed positions inside a launch vehicle fairing. 

The chief limiting factor in accessing space today is the prohibitively large launch cost per unit 
mass of spacecraft. Therefore, the reduction of spacecraft mass has been a primary design 
driver for the last several decades. The traditional approach to the reduction of spacecraft 
mass is the optimization of actuators and structures to use the minimum amount of material 
required for support, deployment, and interconnection. Isogrid panels, composite honeycomb 
panels, and gas-filled inflatable beams all reduce the mass of material necessary to build a truss 
or panel, provide separation between elements, or otherwise apply surface forces to a 
spacecraft structure. An alternative to these "traditional" methods is the use of 
electromagnetic body forces generated by HTSs to reduce the need for material, load-bearing 
support, and standoffs on spacecraft by maintaining spacing, stability, and position of elements 
with respect to one another. 

HTS structures present an opportunity for significant mass savings over traditional options, 
especially in larger systems that require massive structural components. Electromagnetic body 
forces generated by superconducting magnets are used to move and position spacecraft 
elements in lieu of traditional structural components, such as telescoping beams, segmented 
masts and inflatables. HTS structures have less mass per unit characteristic length of the 
spacecraft than aluminum beams and therefore offer the performance benefits of larger 
deployed structures while enabling the stowed structure to fit into existing launch vehicle 
payload fairings. However, the major cost of using HTS structures is the need to cool them to 
low temperatures so that they become superconducting, which requires passive cooling 
structures like heatshields or active cooling subsystems like cryogenic heat pipes. This work will 
also discuss the use of non-superconducting conductors for smaller forces or distances when 
passive cooling is not available. HTSs (which in general are superconducting at temperatures 
below 77K) and room-temperature conductors can be utilized in tandem to perform more 
complex operations. 

This section explains the five primary benefits that HTS structures can offer to the aerospace 
community: 

1. Reduced mass 

2. Larger structures with same launch vehicles 

3. Vibration- and thermally-isolated structures 
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4. Staged deployment, in-space assembly with fewer operations, and partial system 
replacements 

5. Reconfiguration of structures after deployment 

1.1.1. Reduced mass 

HTS structures can change the way mass scales with size; above some characteristic size for 
each spacecraft structure, the HTS system weighs less than other methods. For small, simple 
structures, the mass associated with single or multiple electromagnets and their power systems 
may be greater than that of traditional or other low mass structural options. In those cases, 

HTS structures may not be an advantageous choice (unless they present an additional benefit to 
the system; see benefits 3-5 below). For larger structures, however, the initial mass cost of 
electromagnets is small compared to the mass of a solid, rigid structure. After the payload, the 
structural subsystem is the second largest average fraction of spacecraft dry mass across all 
mission categories, as shown in Table 1, making structural mass the obvious target for mass 
reduction. [1] 


Table 1: Top 4 most massive subsystems average % of dry mass [1] 


Subsystem (% of 
Dry Mass) 

No Prop 
(%) 

LEO Prop 
(%) 

High Earth 
(%) 

Planetary 

(%) 

Payload 

41% 

31% 

32 % 

15% 

Struct/Mech 

20% 

27% 

24% 

25% 

Power 

19% 

21% 

17% 

21% 

Attitude D&C 

8% 

6% 

6% 

6% 


One goal of our work, though not addressed more than qualitatively in the Phase I work, is to 
determine at what size HTS structural options become less massive than alternative structures. 
In order to compare the different alternatives, we define linear structural density as being the 
mass of the structural subsystem (including the mass of those elements of the electrical and 
thermal subsystems dedicated to structural deployment and support) divided by the length of 
the structure (for long, thin structures like booms) and areal structural density as being the 
mass of the structural subsystem divided by the area of the structure (for large, broad 
structures, like heatshields). 

In contrast, the term net areal density is used in this work to refer to the total mass of a 
spacecraft per unit area of a specific characteristic structure, the area of which dictates a key 
performance parameter of the mission, such as the primary mirror of a telescope. Examining 
the net areal density of space telescopes over time reveals an unsurprising trend over time 
towards larger mirrors with lower areal densities. This trend is shown in Figure 1. 
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Figure 1: Chronological progression of net areal densities of space telescopes; mirror diameters listed beside 

points 
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As shown in Figure 1, each successively newer mission has a lower net areal density than the 
last, and each, with the exception of Hubble, has a larger primary mirror than the last. The 
James Webb Space Telescope (JWST) is the first space telescope being built with a primary 
mirror larger than the diameter of its intended launch vehicle and, as such, the mirror must be 
hinged, folded, and then deployed from its stowed position within the launch vehicle. As the 
available launch vehicle fleet remains relatively static in fairing size over the foreseeable future, 
larger space structures will require additional, likely complex and operationally risky 
deployments in order to fit inside the fairing envelope, which suggests another potential 
benefit of HTS structures. 

1.1.2. Larger structures with same launch vehicles 

Many space structures have performance benefits at larger sizes, but spacecraft size is limited 
by fairing envelope dimensions and maximum takeoff weight of launch vehicles. The high 
compaction ratio of HTS structures (related to how the stowed dimensions compare with the 
fully-deployed size of the spacecraft) means that spacecraft designers can reap the benefits of 
larger structures while being less size-limited by the launch vehicle. For instance, a stack of HTS 
coils and tethers that can deploy to a boom length of tens of meters may be less than a meter 
in stowed stack height, whereas aluminum beams can only be as long as the maximum 
dimension of the fairing envelope before requiring hinges and actuators. Figure 2, from 
Bearden, 2001 [2] , shows how launch vehicle costs increase with both orbit altitude and 
satellite dry mass, and how spacecraft going beyond LEO cost twice as much per kilogram as 
those launching to LEO. If HTS structures can reduce the mass of a spacecraft by even 5%, the 
potential savings, both in cost and usable launch vehicle capacity, are enormous. 
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Figure 2: Launch vehicle cost with respect to orbit altitude and satellite dry mass. [2] 
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Example performance benefits of bigger structures: 

• Larger primary telescope mirrors can observe objects farther away because of finer 
angular resolution and increased effective aperture. 

1.22 A 

sinB= — 

The above equation shows the relationship between diffraction-limited angular resolution 
0, wavelength A, and aperture diameter D for a circular aperture. Angular resolution in the 
telescope case is the minimum angle between distinguishable objects in the image 
produced by the telescope's objective, or primary mirror, which serves as the aperture 
when using the above equation. The larger D is, the smaller sin 0 is, which for small angles 
means a smaller 0, thus giving a finer angular resolution and the ability to distinguish 
smaller distances in the telescope's image, unless the payload's resolution is otherwise 
limited by the design of the focal plane sensor. 

• Larger solar sails provide more thrust via greater surface area over which solar pressure 
acts. 

• Larger parabolic radio frequency antennas have higher gain and can enable more distant 
missions or increased transmission data rates. 

AnrjA 


where p is the antenna efficiency, A is the wavelength of the signal being transmitted, and A 
is the effective antenna area. Gain varies directly with A. 

• Larger solar panels can hold more photovoltaic cells and thus generate more power. 

• Longer synthetic aperture radar arrays achieve finest resolution in both cross- and along- 
track dimensions by minimizing along-track (AT) dimension and maximizing cross-track 
(CT) dimension. 

• Larger heatshields can keep more (or bigger) equipment cold. 

While larger structures provide performance benefits, launch vehicle constraints limit 
payload linear dimensions and mass. Table 2 shows the payload capacities of the largest 
currently available and proposed launch vehicles (sans the SLS, which has no hard specifications 
available at this date but is rumored to have an 8m diameter payload fairing), including the 
maximum mass of payload to LEO, the fairing diameter and usable envelope diameter within 
the fairing, and the maximum fairing height for the vehicle. Green rows are the current heavy 
launch fleet; blue are proposed or in-development vehicles; and brown are retired vehicles. 
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Table 2: Current and proposed launch vehicle specifications 


Large Launch 
Vehicles 

Fairing diameter (m) / 
usable envelope (m) 

Max fairing 
height (m) 

Payload to LEO (kg) 

Delta IV Heavy 

5.1/4.57 

9.72 

23000 

Ariane V 

5.4/4.57 

17 

21000 

Proton 

4.35/4 

15.25 

20700 

Atlas VHLV 
(proposed) 

5.4/4.57 

13 

29400 

Falcon Heavy 
(in dev) 

5. 2/4.6 

11.4 

53000 (~26000 if fully 
reusable) 

Space Shuttle 
(retired) 

~4.6 

8 

24400 after Challenger, 
29000 was original 
specification 


Structures larger than about 4.5m in two or more dimensions are thus unable to fit into any of 
the currently available launch vehicles without needing to be condensed in at least one 
dimension. HTS structures provide a means for increasing the deployed-to-stowed length ratio 
of a structure that can be launched by one of the current fleet of launch vehicles. 

1.1.3. Vibration- and thermally-isolated structures 

Structures that are electromagnetically supported (and thus have empty space between them 
instead of solid structures) minimize the conductive pathways between spacecraft elements. 
This enables better thermal and vibrational isolation for sensitive components than if other 
parts of the vehicle are physically joined together. 

The pointing accuracy and image quality of precision instruments is affected by coupled 
vibrations or jitter in the structure from internal or external disturbances (the latter including 
gravity gradient torqueing and solar pressure). Longer structures have a greater displacement 
at their far end than do shorter structures for the same angle rotation as well as a lower 
fundamental frequency, which is undesirable due to potential coupling with the launch vehicle 
during launch and with the attitude control system during operations. Attitude control of the 
spacecraft becomes much more difficult when the fundamental frequency is low. 
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The cryocoolers used with cryogenic heatpipes for cooling superconducting coils generate their 
own disturbance profiles, so vibrational isolation is more approachable in systems already 
possessing passive areal cooling in the form of a sun- or heatshield. 

1.1.4. Staged deployment, in-space assembly, and partial system replacements 

In combination with electromagnetic formation flight, which will be discussed in greater detail 
in Section 2, HTS structures minimize or eliminate connections between spacecraft parts such 
that a heatshield and a mirror assembly of a spacecraft like JWST, for instance, do not have to 
be launched at the same time. Additionally, the entire spacecraft would not have to be replaced 
should one of the two fail because the parts are independent and do not require much or any 
physical assembly. This could give large, complex spacecraft longer lifetimes and the ability to 
be serviced or replaced at distances as great as Earth-Sun L2. 

1.1.5. Reconfiguration of structures after deployment 

Reconfigurability is a spacecraft function that is not currently cost effective since two different 
functional designs on a single vehicle require a large amount of additional 
structural/mechanism mass, and the benefits are usually not worth the added design 
complexity. HTS structures allow for dynamic changes in boom length, solar array placement, 
and heatshield angle, as well as reversible deployments in some configurations - the capability 
to reconfigure could create a whole new paradigm of multi-purpose spacecraft. Much of the 
analysis done in Chapter 4 supports the capabilities of HTS structures for use in 
Reconfigurability, and while HTS structures carry a significant power and thermal burden in 
their use, in large, high-performance systems, unique capabilities like Reconfigurability might 
be worth the trade, especially in designs where the coils do not have to be superconducting and 
in use at all times. This Phase I study and report aims to provide information about the 
potential of HTS structures for such unique functions as well as information about the feasibility 
of HTS structures overall. 
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1.2 Study objectives/Research questions 

The three major questions that this report addresses are related to the technical feasibility of 
HTS structures, the spectrum of viable applications to spacecraft, and the new operations 
enabled or made substantially more realistic by HTS structures respectively. 

Question 1: Can we use electromagnetic forces generated by and acting between high- 
temperature superconductor current-carrying coils to move , unfold, and support parts of a 
spacecraft from its stowed position? 

In order to answer that question, we first define functional configurations of coils that can be 
used to perform specific operations onboard a spacecraft, distinguished from each other by 
their degrees of freedom, starting and ending positions, number of coils, and whether they are 
done at initial deployment of a structure or at a point later in the operational lifetime. 

HTS coils can repel, attract, and even shear with respect to one another; in a flexible, non- 
circular coil, elements of the same coil can perform these actions on one another, deforming 
the shape of the coil over time. In combination with various boundary conditions, these 
operations lead to the seven functional coil configurations, depicted in Table 3, that we have 
identified for use in spacecraft deployment and support activities: four (1-4, in blue) are for 
initial deployment and three (5-7, in green) are variants of the deployment configurations for 
use in overall spacecraft shape change during the spacecraft operational lifetime. Table 3 
describes the configurations and the spacecraft operations they could perform. These 
configurations will be capitalized when referred to in this report. 
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Table 3: Seven potential functional coil configurations 


Configuration 



Description 

Uses 

A single, flexible HTS coil that is 
folded in its stowed position 
and uses its own magnetic field 
when current is run through it 
to Expand to flat shape 

To deploy and hold taut the 
perimeters of large 
membranous or flexible 
structures 


2. Inflate 


Deployment 



A 3D structure is built with two 
or more repelling coils in a 
configuration that creates a 
space between the two, 
Inflating a structure (bounded 
by flexible walls or tethers) 


To create a volume (such as 
a tank or toroidal perimeter) 
or to increase surface area 
(for solar cells) 


3. Unfold 

Deployment 

A series of coils embedded in or 

To deploy long, flat 



attached to a structure that is 

structures (like solar panels) 



stowed folded and must be 

or to take advantage of 



Unfolded to become 

mechanical resistance at 



operational (folds can be 

folds to create variable 



hinges, springs, or couplings) 

angles 


4. Separate 


Deployment 



Corollary of Inflate, in which 
two or more coils repel each 
other in series facing each other 
to Separate two parts of a 
spacecraft; tethered or 
membranous structure 
connecting 


To put large, controllable 
distance between two 
sensitive parts of spacecraft 
(such as a nuclear reactor, 
astronauts, optics, 
thrusters) 


5. Deform 


Operational 


* °o° 

o * ° 


Two or more coils embedded in 
parts of the structure act 
magnetically on each other to 
temporarily Deform or change 
the shape of the spacecraft 


To reduce radar cross 
section (RCS) or adjust 
shape for avoidance of 
debris 


6. Reconfigure 


Operational 



Corollary of Deform, except 
Reconfigured state is 
sustainable and lasting 


To reduce drag profile or 
Reconfigure satellite for 
different ConOps 


7. Refocus 


Operational 



Two or more coils adjust their 
magnetic state such that an 
antenna or mirror is Deformed 
to Refocus it. Orientation of 
coils dependent on original 
shape of mirror or antenna. 


To change focus lengths of 
mirrors and gains of 
antennas by reforming or 
moving their dishes, 
mirrors, or horns to Refocus 
them on a new target 
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The four deployment configurations (1-4 in Table 3) are numerically modeled in MATLAB and 
Simulink to verify that the deployments can be performed as described in Table 3 and are 
validated in the far-field (where the coils are very far apart from one another) against the 
dipole approximation of coil forces and torques. Chapter 3 describes these models and how 
they were constructed and validated. 

Question 2: For which operations does this technology represent an improvement over existing 
or in-development options? 

In order to investigate the performance of HTS structures versus other low-mass structural 
technologies, the ranges of key performance parameters over which HTS structures are viable 
options for spacecraft design must first be determined. Chapter 4 characterizes the 
performance of HTS structures using the models introduced in Chapter 3, while in Chapter 5 
they are compared to other structural options, both traditional rigid metal structures and 
"alternative" structures made with different materials and deployment operations. 

Question 3: What new mission capabilities does this technology enable? 

The Reconfigure, Deform, and Refocus operational configurations all represent mission 
capabilities that have heretofore required structures and mechanisms that are too expensive in 
metrics like mass, size, or power or are too complex and therefore too risky to justify for the 
additional performance benefits that they offer. This is why spacecraft do not commonly 
include these capabilities. 

Using magnetic forces, however, such shape-changing mission capabilities are not structurally 
much different from the deployment and support configurations required to move from a 
stowed configuration to a deployed one, so including the ability to Reconfigure or Deform 
utilizes the existing structural architecture more efficiently than in non-electromagnetic 
structures. Operational configurations will be discussed further in Chapters 3 and 4. 

The goal of this report is to provide analyses of the dynamics of rigid (and, briefly, flexible) 
electromagnetic coils when exerting forces and torques on one another, such that we can 
support a positive answer for Question 1, narrow down the field of potential deployment and 
support functions that are worth further study for Question 2, and explore new spacecraft 
architectures that have not previously been worth seriously entertaining for Question 3. 
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2 Chapter 2 - Background 

This chapter introduces the foundational scientific principles for HTS and reviews previous 
development efforts to mature key enabling technologies for HTS electromagnetic structures. 
We discuss the enabling scientific principles and phenomena that form the basis of our study, 
draw lines between HTS structures and previous space electromagnetic work, and briefly 
examine the maturity of the HTS structure concept before moving into more theoretical 
analysis of the deployment dynamics in Chapter 3. 


2.1 Scientific principles enabling HTS structures 

HTS structures are enabled by the fundamental scientific principles of electromagnetism and 
the unique environment of space, as well as the developments made over the last several 
decades in manufacturing superconductors and controlling electromagnets onboard multi- 
system vehicles. The industry development is explored in later sections; this section describes 
the underlying fundamental scientific principles that enable HTS structures, including: 

1. The creation of Lorentz and Laplace forces via interaction of a magnetic field and 
current 

2. The Meissner effect, superconductivity, and the manufacturing of HTS wire 

3. Enabling characteristics of space environment (microgravity and vacuum) 

2.1.1 Generation of Lorentz and Laplace forces 

HTS structures operate using electromagnetic forces to push, pull, and move with 
respect to each other. Electromagnetic forces (called Lorentz or Laplace forces depending on 
whether the force is acting on a single charge or a current of charges respectively) result from 
the interaction of a magnetic field and a current. 

A point charge q moving with velocity v in external magnetic and electric fields, B and E 

— > 

respectively, experiences a Lorentz force F, given by: 


F — q[E + v x B] 


where qE is the electric force, and qv x B is the magnetic force. The macroscopic force on a 

wire is the magnetic force and is called the Laplace force. It is generated by a magnetic field 

— > — > 

B on a wire carrying current / (a stream of point charges) as follows: 

F — I I dlxB 

wire 
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The magnetic force is thus orthogonal to the wire and to the orientation of the magnetic field at 
the point of calculation. In order to calculate the magnetic field B for use in determining the 
Laplace force on a wire or coil of wire, one can use the Biot-Savart law, which can be derived 
from Ampere's law and Gauss's law, to compute the resultant magnetic field vector B at a 
position r with respect to a steady current /: 



Idl x r 


|r| 


3 


wire 


In the deployment modeling work that will be described in Chapter 3, the magnetic field and 
the Laplace forces across a current-carrying wire are approximated over time by implementing 
the Biot-Savart law numerically, discretizing electric current elements in order to determine the 
magnetic field generated by arbitrary configurations of rigid (meaning a fixed, non-changing 
shape) and flexible coils. Knowledge of the magnetic field at each point in space around a 
current-carrying wire allows the calculation of the resultant force upon another current- 
carrying wire as a result of that magnetic field, which can then be used to determine the 
number of turns a coil requires or how much current it needs to carry in order to deploy a 
structure. 


2.1.2 Meissner effect, superconductors, and manufacturing of HTS wire 

Superconductors enable HTS structures because they are able to generate much larger forces 
via their larger current carrying capacity, which increases the distance over which they can 
work for the same amount of mass. Superconductors are materials that conduct electrical 
current perfectly below a critical temperature T c ; superconductors have zero resistivity, with 
negligible quantities right around their T c . Any current through them will persist significantly 
longer than through a non-superconductive material. Superconductivity is characterized by the 
Meissner effect, the expulsion of an external magnetic field from a superconductor once cooled 
below its T c during its transition to a superconducting state. Every superconductor has a critical 
temperature, external magnetic field strength, and current density above which 
superconductivity ceases, shown graphically in Figure 3. 
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Figure 3: Critical surface for Type II superconductor [3] 

High-temperature superconductors, or HTSs, are those superconductors with T c s above 77K, or 
the boiling point of liquid nitrogen, enabling them to be cooled to a superconducting state 
using LN2. The higher the T c , the less it costs (in terms of power, storage, and consumables, for 
applications where a cryogen is not recycled) to do the cooling. HTS development and the 
subsequent development of HTS wire has led to a broad array of applications for 
superconductors, including long distance power transfer, electromagnets, and energy storage. 

There are two types of superconductors; Type-1 only exhibit the Meissner effect with one 
critical field strength above which superconductivity ceases. Type-ll superconductors, which 
include all high-temperature superconductors, as well as some low-temperature 
superconductors (LTSs) with T c s too low to qualify as "high-temperature", also exhibit a 
"mixed-state" Meissner effect that increases their critical magnetic fields and configuration 
stability. Because of this effect, Type-ll superconductors are often used in superconducting 
magnets in the form of coils of wire made with superconductor filaments embedded in support 
material less than a millimeter in diameter. 

The "mixed" Meissner effect is different from the Meissner effect in that some magnetic field 
penetrates the superconductor through filaments of normal-state material, and the material 
can support higher magnetic fields before superconductivity breaks down. There are thus two 
critical field strengths in Type II superconductors: beyond the first field strength, where 
superconductivity would cease completely in a Type I superconductor, a vortex ("mixed") state 
exists in which some magnetic flux is allowed to penetrate the material while it continues 
superconducting. Beyond the second, higher critical field strength, superconductivity ceases. 

Type-ll superconductors are the only type of superconductor used in wire. Many of these 
types of wire are made with HTSs to lower cooling costs, especially for non-magnetic 
applications. Some wire, however, especially that which is used in powerful electromagnets 
like those in the Large Hadron Collider, is made with LTSs that need to be kept much colder but 
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can sustain much higher current densities than HTSs. For example, niobium-tin (Nb3Sn) has a 
T c of 18. 3K and can withstand magnetic field strengths up to 30 tesla (with a record current 
density of 2643 A/mm A 2 at 12T and 4.2K) [4], HTS wires do not have the same current 
densities, and, as a result, they cannot generate magnetic fields as high as LTSs. But, though 
pure HTS materials are just as brittle as LTS materials, HTSs can be constructed into flexible, 
durable wire. The flexibility of HTS wire enables the Expand configuration that we proposed 
earlier in this section for deployment of a single, folded and stowed coil into a large, flat, 
expanded coil. 

Compared to standard, room-temperature conductors, HTS wires are able to create larger 
magnetic fields and sustain higher current densities, with little-to-no resistive losses through 
the wire (compared to high resistive losses in copper and aluminum). While room temperature 
conductor coils can be used to magnetically operate on each other, the Laplace forces able to 
be generated on each other are significantly lower than those that can be created with HTS 
coils, due to the lower induced magnetic field. HTS wire thus enables electromagnetic 
structures with multi-meter separation between coils, which in turn enables larger vehicles and 
the performance benefits enumerated in Chapter 1 that accompany larger structures. 

Whereas several decades ago it was difficult for manufacturers to create HTS wire of any useful 
length, HTS can now be created in lengths upwards of 1500m [5], reducing the need for splices 
that result in resistive losses of power to heat. Individual strands of tape-shaped (wide and 
thin) cable are multilayered with various substrates and insulators for mechanical strength and 
chemical stability while at the same time trying to maximize the current density of the wire. 
Figure 4 shows a cutaway view of an HTS wire. 


Copper Stabilizer 


20 pm 


Silver Overlayer 

(RE)BCO - HTS (epitaxial) 



Figure 4: Cutaway view of a 2G HTS wire [6] 


Roebel cables, shown in Figure 5, implement a way of winding together individual wire strands 
to reduce AC losses due to self-field interaction with the current running through each strand, 
reducing the critical current I c able to pass through that wire. When determining the 
configuration of turns in a coil or strands in a cable, it is especially important to consider how 
the self-field affects the I c in each strand to optimize for the maximum current density in the 
coil cross-section possible. 
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Figure 5: 15-strand HTS Roebel cable [7] 

Superpower Inc, a subsidiary of Furukawa Electric Co., Ltd, manufactures a flexible "2G", or 
second generation, HTS wire with a minimum bend diameter of 11mm at room temperature 
(and at least 25.4 mm at superconducting temperatures) and maximum hoop stress of 700 MPa 
before irreversible degradation of the critical current. Recent axial stress-strain measurements 
conducted at the Naval Research Laboratory on standard 4 mm wide SCS4050 production wire 
gave a yield stress for the wire of 970 MPa at 0.92% strain. [8] The 2G HTS wire can carry 
between 250 and 350 A/cm wire width (or 25 - 35 kA/cm2 with a ,1mm wire thickness) through 
650m lengths (80-110A critical current with a 4mm wire width) [6], We use the properties of 
this 2G wire in our models. 

2.1.3 Space environment 

On the ground, gravity and the need to cool superconductors in an otherwise room- 
temperature environment makes using electromagnets as actuators or structural support 
difficult and not broadly useful. In the space environment, however, an electromagnet does not 
need an enormous magnetic field to actuate components, and there is no air transferring heat 
into the magnet by convection, making space a potentially favorable environment for the use of 
HTS structures. 

The microgravity environment of space enables HTS structures because spacecraft elements 
can be actuated without overcoming gravitational forces. Thus, only small forces are needed to 
cause motion or actuation, reducing the necessary size and current of electromagnets used for 
such tasks compared to what would be needed, for instance, to repel a coil upwards on the 
Earth's surface. 
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The vacuum of space is both beneficial and detrimental to the use of superconductors. If it can 
avoid radiated input from the sun. Earth, or the space vehicle, a thermally isolated HTS 
subsystem can maintain superconducting temperatures without cryogens. However, it is 
difficult to remove heat from the system in a vacuum if it cannot be isolated from conductive or 
radiating heat sources. This thermal environment makes more advanced cooling systems, like 
the cryogenic heatpipe described in Section 2.2, necessary for maintenance of the 
superconductor below its critical transition temperature. 

2.2 Enabling technology and previous work 

A number of previous and ongoing industries and projects lend credibility to the concept of 
using HTS to actuate spacecraft elements with respect to one another, including (1) dynamics 
and control technologies and testbeds like MIT SSL's electromagnetic formation flight (EMFF) 
testbed and the University of Maryland's (UMD) RINGS project (slated to fly to ISS in December 
2012) Peck's flux pinning [9], and Pedreiro's Disturbance Free Payload [10]; and (2) thermal 
control innovations like Sedwick's cryogenic heatpipe for HTS coil cooling. 

1. Dynamics and control of a spacecraft without a support structure become a significant 
challenge and represent the majority of the work that must be done to design a system of 
electromagnets, since actuation and support are dictated by the feedback-informed 
distribution and direction of power through wire. The MIT SSL and UMD Space Power and 
Propulsion Laboratory have studied the problem of controlling electromagnetic coils in 3- and 
6-DOF systems for the last decade, showing that control of free-flying or tethered vehicles with 
electromagnets is feasible and reinforcing the concept of using HTS coils to deploy a structure. 

Electromagnetic formation flight (EMFF), has been the subject of study by the MIT SSL since 
2002, funded by the original NIAC program (Phase I and Phase II studies), NRO Dll, the JPL TPF 
Program, NASA GSFC SBIR, and DARPA. In 2003, the MIT AeroAstro senior capstone design 
class focused on the development of an algorithm for formation control using the inter-vehicle 
forces from steerable electromagnetic dipoles. As part of the class, and with funding from the 
NRO Dll and the NASA Institute for Advanced Concepts, the MIT group implemented a 
laboratory testbed that used high temperature superconducting (HTS) wire to create 
electromagnetic coils. In a follow-on effort, MIT and Aurora Flight Sciences (AFS) investigated 
the performance of a room-temperature conductor version of EMFF (termed microEMFF). This 
effort included both a performance assessment and construction of a ground-based microEMFF 
testbed. Resonant Inductive Near-field Generation System (RINGS), is a joint UMD SPPL/MIT 
SSL/AFS project that, among other things, is a microEMFF demonstration onboard the 
International Space Station. These programs have included intensive control algorithm 
development as well as TRL development of the HTS EMFF concept to TRL 4 and the microEMFF 
concept to TRL 5, once RINGS begins operations on ISS. 

The basic concept of EMFF is to provide actuation in relative degrees of freedom for formation 
flight systems using electromagnetic forces/torques and reaction wheels. The positions of 
system elements can change so long as the center of mass of such a system remains fixed 
(barring external input). The motivation for EMFF in multi-vehicle systems is fourfold: 
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• Station-keeping for distributed satellite systems 

• Replacement of consumables (thrusters) 

• Elimination of thruster plumes 

• Enabling of high delta-V formation flying missions 




Magnet and 
cryogenic 
containment 


Metrology 


Reaction 

wheel 


supply 


Power 
Boards & 
Batteries 


EMFF Testbed Parameters 

Mass (single vehicle) 

36 kg 

Height 

1 18m 

Maximum Diameter 

0 90 m 

Operating duration 

- 20 minutes 

Table dimensions 

3 m x 15 m 


Figure 6: MIT SSL EMFF Testbed [11] 


In the 3DOF EMFF testbed built at MIT SSL, shown with subsystem labels and with both test 
vehicles in Figure 6, vehicles exert electromagnetic forces and torques on each other using two 
perpendicular FITS coil electromagnets each. A 6DOF system would use three orthogonal HTS 
coils to create a completely electromagnetically steerable magnetic dipole. The reaction wheel 
decouples torques, allowing independent control of a vehicle's rotation. 


During MIT SSL's previous work with this experimental testbed, several control algorithms were 
developed and validated. These results show that the relative position of spacecraft can be 
controlled and stabilized using 1 to 3 electromagnetic coils and 3 reaction wheels per vehicle 
[12]. While this particular testbed was not a consumable-free, closed-loop system (LN 2 cooled 
the superconductors but was allowed to boil off), the cryogenic heatpipe previously discussed 
enabled a closed-loop system with no LN2 loss. 


RINGS, or Resonant Inductive Near-field Generation System, is an electromagnetic formation 
flight and wireless power/data transfer testbed that will begin operations on the ISS in 
December 2012. It is an add-on to the SPHERES (Synchronized Position Hold Engage and 
Reorient Experimental Satellites) developed by MIT SSL and is currently being operated by 
NASA Ames Research Center as a National Lab onboard the ISS. RINGS is an advancement of 
the microEMFF concept in that it uses room-temperature conductors instead of 
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superconductors, and its use on the ISS will advance the TRL of microEMFF to 5 via brassboard 
demonstration in a relevant environment. The RINGS consist of two 155-turn, 0.685m outside 
diameter coils of non-insulated air-cooled aluminum cable enclosed in a protective plastic 
shielding. One coil is mounted on each of two SPFIERES as shown in Figure 7 and, using the 
SPFIERES' C02 thrusters to steer the dipole, the coils can be used to maneuver with respect to 
one another magnetically. 



Figure 7: A SPHERE outfitted with RINGS hardware in laboratory 

The primary differences between the application of FITSs in EMFF and in structures/deployment 
are (1) the presence or absence of physical connection paths between interacting coils and (2) 
the level of real-time direction the system has over the movement of coils. 

EMFF vehicles can steer their dipoles in response to real-time system input to maneuver with 
respect to one another and have reaction wheels to cancel out torque when shearing with 
another vehicle. FITS structures are limited by their physical connection to the elements that 
they are moving, and their movements are planned for and constrained to certain, predictable 
paths, as is desired when deploying spacecraft elements to operational configurations. It is 
important to note that FITS structures can isolate different parts of a spacecraft vibrationally to 
a greater extent than can solid structures, but unless the elements are formation flying, there 
are still some transmission paths of forces and vibration, limiting movement of a coil in at least 
one degree of freedom. The EMFF project did not focus on deployment dynamics and the 
stability of transient states from one position to another without capacity for 6DOF control and 
steering of the resultant dipoles, thus further work is required to characterize this feasibility risk 
(and control design challenge) for the FITS structures application. Although controllability is a 
major component of technical feasibility for FITS structures, it is not the focus of this report. 


Electromagnetic actuation on spacecraft has been the subject of research by a number of 
groups because of its potential for vibrational isolation and reconfiguration. Disturbance Free 
Payload, or DFP, is an architecture developed by Pedreiro et al for use in vibration isolation of 
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spacecraft payloads from the rest of the spacecraft using multiple electromagnetic actuators 

[ 10 ]. 

Flux pinning is an electromagnetic phenomenon in superconductors that enables highly 
accurate positioning between an electromagnet and a permanent magnet or other magnetic 
field source that can be used for safer and easier in-space assembly and reconfiguration of 
space systems in orbit. Small defects in the superconducting material can increase critical 
magnetic fields and contribute to the stability of the system by fixing vortex points, or 
filaments, through which magnetic field lines pass and are subsequently pinned when the 
superconductor is cooled below its T c . Pinning of the field lines pins the source of the magnetic 
field in position and orientation. A number of joints, hinges, and interfaces can be created 
using configurations of permanent magnets to limit or direct the motion of a superconducting 
cube about the permanent magnets. As of 2009, multi-centimeter gaps between the magnets 
were supportable with high stiffness and damping with small masses. Flux pinning as a means 
of supporting and manipulating space structures is the subject of work by Peck et al at Cornell 
University [9], 

Stiffening of membranous or low net areal density structures has been proposed on multiple 
occasions, such as Zubrin's "magsail" concept in which a large membranous plasma wind sail is 
held taut using a flexible HTS coiled around its perimeter to repel itself into a circle and keep 
the sail under tension. Zubrin ultimately selected a non-magnetic deployment system (rotating 
booms to deploy the sail initially using centrifugal force), citing "reliable deployment" as a key 
issue for magsails. [13] A previous NIAC study performed by Powell et al [14] proposed a system 
like that being studied in this report: flexible cables made of high-temperature 
superconducting wire Expanding or Inflating in order to serve as deployment actuators, 
perimeter support, and standoff structures simultaneously for large scale spacecraft. Powell 
primarily discusses the HTS wire properties and materials and does not go into detail on the 
deployment process. This report will explore such magnetic deployment in detail in an effort to 
determine if magnetic deployment is feasible, and if so, for what types and sizes of structures. 

Other past studies of deployment or support using electromagnetism in membranous 
structures include microwave beam-driven spin deployment of solar sails (Benford, [15]) and 
membranous structures containing conductive meshes that can be shaped using magnetic 
pressure from permanent magnets or electromagnets (Amboss, [16]). 

For applications that require long, low-mass structures, HTS structures are more suited than 
flux pinning or DMP to deployment and support in those situations. The possibility remains, 
however, of using a combination of flux pinning, DMP, EMFF, and HTS structures to 
magnetically assemble, deploy, and support large and complex space structures. Other light 
structural technologies like inflatables or tensegrity structures do not have active control over 
their shape; they attempt to deploy and either succeed or fail, and then they passively maintain 
their configuration for the lifetime of the vehicle. The ability to control and change the 
magnetic field means that electromagnetic structures are able to change their shape after 
deploying, depending on their boundary conditions. 
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2. Thermal control down to cryogenic temperatures is necessary to induce 
superconductivity in HTSs. This level of thermal control is not normally present on spacecraft 
unless the vehicle has a payload that requires special cooling, making the need for cooling one 
of the primary disadvantages of HTS structures. Additionally, HTS performance is dictated by its 
warmest temperature, meaning that cooling devices are required at every point along the coil, 
which is not feasible (or even accomplishable), for large coils. Cooling equipment that can cool 
the full perimeter of a coil is needed and must be scalable to different size coils. Such 
equipment should preferably be consumable-free so as not to limit the lifetime of the mission. 

If the thermal control requires too much mass and power, the thermal subsystem could make 
HTS structures less competitive with traditional structures by nullifying any potential savings of 
HTS structures. 

The cryogenic heatpipe developed by Sedwick and Kwon [11] for use in EMFF, pictured in 
Figure 8, provides a starting point for continued development and improvement of a 
competitive HTS coil cooling system by accomplishing consumable-free isothermalization to 
cryogenic temperatures (the same temperature at every point on the coil) with a single 
cryocooler. 



Figure 8: Cross-section of heat pipe (left) and cryogenic heatpipe (copper casing) in open toroidal vacuum 

chamber (right) 


This rigid heatpipe uses nitrogen as a working fluid and a stainless steel mesh as the wicking 
structure, with the HTS coil enclosed within the heatpipe instead of residing externally as with 
traditional heatpipes. A condenser, operated with an LN2 reservoir (like in the laboratory setup 
shown in Figure 8) or powered by a cryocooler, extracts heat from the system while the mesh 
passively wicks the condensed LN2 around the pipe via capillary action on several layers of wire 
mesh to cool and achieve isothermalization of the coil. Isothermalization means that minimal 
temperature gradients exist across the entire coil. Isothermal conditions are desired because 
the warmest point in the coil is the limiting factor in the performance of the HTS wire. When 
used with a cryocooler, as in the composite heatpipe design discussed next, this heatpipe is a 
closed system wherein the working fluid is conserved so long as the system remains sealed. 
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UMD SPPL conducted additional heatpipe development under a NASA STTR, in cooperation 
with Axis Engineering Technologies, Inc. [17]. The principal heatpipe design remained the same; 
however, the focus was on manufacturability of the heatpipe, which was previously found to be 
problematic. The issue was that the superconducting coil had to be a continuous spiral, 
requiring that the housing be formed around the coil. The housing must also be thermally 
conductive but ideally electrically insulating, since eddy currents developed in the housing 
would generate heat and dissipate power. 

The resulting design is shown in Figure 9. The clamshell design is constructed from a thin- 
walled composite material, instead of the copper used in the first heatpipe, held together with 
a set of inner and outer clamp rings. The composite material is electrically insulating, but the 
thin wall allows for good heat conducting radially. The outside of the pipe is subsequently 
coated with MU. 


The cryogenic heatpipe, about two meters in diameter with 100 turns of HTS coil within its 4cm 
pipe diameter, was designed with a heat capacity of 100W (determined from the measured 
radiated heat input from the vacuum chamber walls) and succeeded in keeping the HTS coil 
below its critical temperature of 110K. The mass of the 2m diameter copper pipe was 
approximately 70.5 grams. In LEO, however, the worst-case heat capacity of a 1 meter 
diameter MU-coated coil (the size used in the EMFF testbed) is 5W [11], 5W of heat can be 
extracted using a relatively low-mass (2.1 kg) cryocooler like the Sunpower CryoTel MT [18]. 


Outer Ri 



Sandwich P*' es 


Figure 9: Two-part composite heatpipe showing an HTS coil inside 
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2.3 Technology readiness 

HTS structures are a logical evolution of the previously realized or studied applications of 
electromagnets in space, and they bridge the gap between electromagnetic actuation, support, 
and positioning as a system that is capable of all three functions. In the future, spacecraft 
designers can tie together many of these positioning and actuation applications into a single or 
collection of multiple spacecraft, utilizing formation flying, easy in-situ assembly, isolation and 
high stowed-to-deployed size ratio structures to enable spacecraft that have never before been 
possible. Therefore, in this report, we will address several of the questions relevant to the 
feasibility and viability of implementing HTS structures through modeling and trade analysis. 
Furthermore, a goal of this work is to mature the technology readiness level (TRL) of HTS 
structures from conceptualization and early analysis of this application (TRL 2) to a more 
detailed analysis and hardware proof-of-concepts (TRL 3). Phase I encompasses the transition 
from TRL 2 to early 3 while our Phase II work aims to mature the technology further through 
TRL 3. 
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3 Chapter 3 - Theoretical Approach 

This chapter will detail the equations of motion that describe the dynamical evolution of 
current-carrying coils in space and how this motion was simulated using MATLAB and Simulink 
models. In this report, we will be primarily discussing the motion of two rigid coils under a 
variety of boundary constraints; we will briefly discuss the motion of flexible coils in this 
chapter with a more in-depth discussion and analysis reserved for our Phase II report. These 
particular models were selected because between the two of them, all seven functional 
configurations presented in Chapter 1 could be modeled given appropriate boundary 
conditions. 

In general, any current-carrying coil will experience forces and torques due to other current- 
carrying coils and magnetic fields in the vicinity. One of the purposes of this work is to 
determine how the mutual evolution of a network of coils is dependent upon the initial 
positions and orientations of the coils and how that evolution can be controlled for deployment 
and subsequent support purposes by imposing various constraints upon them, such as flexible 
tethers and hinges. Accurate models of coil dynamics can be used to characterize the range of 
system parameters over which electromagnetic structures are potentially useful. In this 
chapter, these models and the methods used in their construction will be introduced for use in 
studying the ranges over which each of the functional configurations is feasible. 

3.1 Dynamics of unrestrained coils 

The coil dynamics can be described fully by Newton's laws of motion. Newton's 2 nd law will be 
employed in both its rectilinear form, for modeling the motion of the center of mass of the coil, 
and in its angular form, in order to describe the rigid body rotation of the coils through space. 
Newton's laws are first supplemented with the appropriate electromagnetic laws introduced in 
Chapter 2 to describe the electromagnetic forces and torques exerted on each coil by one 
another in a six degree-of-freedom model and later with mechanical forces that reflect the 
influence of constraining devices (tethers and hinges) on the coils. 

3.1.1 Rectilinear motion 

Consider a set of current loops, one of which we label by the subscript /, where / = 1,2,..., N, 
with N being the number of coils in our system. The magnitude of the current in coil / is given by 
/j. At some point P in space (defined with respect to a suitably-chosen inertial coordinate 
system), the magnetic field due to a directed, differential coil element df is given by Ampere's 
law: 



pli dlj x f Pi 
An |r Pi | 3 


where f Pi = rj — f P is the vector from P to the differential coil element, dl t is the vector 
between two infinitesimally close points on the current loop, and p is the magnetic constant, 
equal to 1.26 x 10~ 6 T-m /A. 
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The total magnetic field at P due to coil i is given by the integral over the entire loop, or: 



f pf dl t x f Pi 
J 47 r \r Pi \ 3 

Loop i 


Now, suppose P happens to be on another current loop (denoted “j"), with the magnitude of its 

— ) 

current given by Ij. Let the differential coil element at P be given by dlj-, then, by the Lorentz 
law, the force exerted on this differential coil element by the magnetic field due to coil / is given 
by: 


dFj — Ijdlj x Bj 


— Ij d lj 



Loop i 


Ph dli x Vji 

4n \y | 3 
I 'ji I 


where r i; - has been rewritten from f iP to emphasize that the magnetic field from coil / is acting 

— > 

to create a force on the differential coil element dlj. I t and Ij are constant values around the 
current loop, and in cases in which each coil is made up of wire turned multiple times, if is 
the number of turns in the /th coil, we may define an effective current I e ff,i as follows: 


hff.i 


= n t I t 


The total force exerted on coil j by coil / is then given by integrating the Lorentz force over coil j: 


F j = 


P^eff.i leff.j 

An 


dlj x 


dh x Vji 


Loop j 


Loop i iJi\ 


By Newton's 2 nd and 3 rd laws, the force acting on coil / from coil j is given by 


Fi Fj 77tjClf, 

these forces acting on the centers of mass of each coil. For a system with more than two coils, 
the force is equal to the summation of the contributions from each of the N- 1 other coils (A/-1 
instead of N because the self-contribution of force from coil j is 0). 


Since a coil responding to the magnetic fields of other coils will move and rotate, the vector 
variables of current element direction and position in the above equation change with time for 
each coil in the system. As such, the delay in the response of a given coil due to the finite 
propagation speed of electromagnetic waves must be taken into account. In our case, however, 
so long as we are interested in time scales much greater than |r i; - |/c, where c is the speed of 
light (3xl0 8 m/s), the effect of the delay will be negligible. 
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By denoting the position of the center of mass of the y'th coil in the inertial coordinate system 
by an overbar, fj, we have: 


-» d 2 rj 

F i= m iMT 


N 

Z P-leff.i hff.j 
An 

7 = 1 
j*l 


dlj x 


Loop j 


dli X T~ji 


Loop i I r n\ 


There is no clean analytical solution for this problem besides the above double integral. For an 
axially symmetric coil, however, the magnetic field at a point P along the axis of symmetry out 
of the plane of the coil has the analytical solution: 



pnlR 2 

2 (R 2 + x 2 ) 3 /2 


where B is calculated in Tesla, R is the coil radius in meters, / is the current in amps, n is the 
number of turns of wire in the coil, and x is the distance in meters along the axis from the 
center of the coil. This solution is useful for calculating the force on a single particle along the 
coil's axis, but not for the force at any other point. Therefore, a numerical approximation of the 
force between the coils is required in order to study the behavior of the coils over time, 
accomplished by numerically performing the line integrals above. In this report, we investigate 
the dynamics of the coils using numerical models constructed in MATLAB and Simulink. 


3.1.2 Rigid body dynamics and rotation 

In addition to the translational, rectilinear motion of two coils in space, the torque and 
resultant rotation induced by asymmetric forces on differential elements across the coil must 
be taken into account. 

By the force expressions given above, the torque exerted on the /th coil about its center of mass 
is given by: 


N 




7 = 1 
j*i 


4n 


( n - rf) x dli X 

Loop i Loop j 


d lj x r i j 


1 1 j I 


where r £ is the location of the differential element dl t with respect to the global inertial 
coordinate system, and (r £ — ff) is the vector from differential element dl t to the center of the 
coil, providing the moment arm of the torque formulation. 
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As with any analysis of rigid body motion, care must be taken in how the orientation of each 
body (in this case, the individual coils) is parameterized. In this report, the orientation of the 
coils in space is parameterized by means of a quaternion. 

The use of quaternions to describe the orientation of a body in space is well established— for a 
full account of the theory, the reader is referred to, e.g. [19]. We describe the approach in 
relation to coil dynamics briefly: 

We attach a coordinate system (called the body coordinate system and denoted (X^j) to the 
ith coil. The orientation of this coordinate system is encoded within the quaternion q t via the 
direction cosines matrix A ir which specifies the coordinates of {X^} with respect to a suitably 
chosen global inertial coordinate system. 


The matrix A ir when multiplied to the right by the set of vectors that define the positions of the 
points around the coil around its center, effects the appropriate rotation upon the coil. The 
dynamical response of the coils is calculated using Newton's 2 nd law in angular form, or 


Ti = [I]idi 


Here, [7]j is the moment of inertia matrix of the ith coil calculated with respect to its center of 
mass, given by: 


U]i = 


yniRi 

2 

0 

. 0 


0 

rriiRi 

2 

0 


0 

0 

rritRi - 


We have chosen so that it is aligned with the principal axes of the coil; the third entry in 
the diagonal of [7]j corresponds with the axis normal to the plane of the coil. The angular 
acceleration of the body is <5j. From this, the evolution of the quaternion— and the orientation 
of the coil— can be advanced accordingly. [19] 


3.2 Incorporation of constraining elements 

The equations of motion presented above are for coils in empty space. Ultimately, mechanical 
elements (possibly in conjunction with appropriate current controllers for the coils) must be 
employed in order to execute desired geometrical configurations of the coils in space. The coils 
and these constraining elements will then serve as the primary structure of the spacecraft. 
Below, we describe incorporation of such constraints into the equations of motion. 

3.2.1 Tethers 
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For the Separate functional configuration, the use of tethers to connect magnetic coils to each 
other and restrict their motion eliminates the need to control currents in order to keep coils 
positioned the desired distances from one another, as they would have to be in 
electromagnetic formation flight. Tethers also reduce the risk associated with two or more 
coils accelerating away from each other during deployment. This risk reduction is performed 
via the tethers providing a corrective "spring" force for any overshoot of the desired separation 
distance. In this report, three tethers equally spaced around the perimeter of the coils are used 
to restrict the motion of the coils, as shown in Figure 10. In future work, we will consider more 
complex tether configurations. 



Figure 10: Diagram of 3-tether spacing between two coils 


Three tethers are used instead of two because two tethers do not restrict the rotation of the 
coil around the axis drawn through the tether attachments (visualized in Figure 11), which 
could lead to a 180 degree flip of one or both coils once the tethers reach their full length 
during deployment. 


\ 


Figure 11: Diagram of 2-tether spacing and unrestricted motion 
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In the Separate deployment, the coils start close to one another with tethers slack, as shown in 
Figure 12. Suppose the coils are aligned along their central axes, and the tethers— all identical 
to one another— are parallel to these axes. Once current is run through the coils in opposing 
directions, the coils begin to accelerate away from each other until the coils are L t apart, where 
L t is the unstretched length of a tether. If the tethers are stretched to a length L > L t , an 
elastic force due to the tension of the tethers is effected. If we define A L = L — L t , then the 
elastic force of a single tether on a given coil is \P e i ast itether \ — kAL, where k is the spring 
constant of the tether with the force directed towards the other coil. If there are three equal 
length and material tethers, then the total elastic force acting on a coil is |F eiast | = 3kAL. In 
this model, the tethers are not rigid and do not provide compressive forces for distances 
shorter than L t . 


L 

Jl 



Thus, in a two coil system with three identical tethers, when L < L t , the forces acting on the 
coils are solely to due to electromagnetic interaction, while, when L > L t , the forces acting on 
the coils will be due to both electromagnetic interactions and elasticity, the magnitude of the 
elastic force given by 3kAL. 

In the more general case of coils and tethers not aligned, the effect of the tethers on the coils 
may be determined as follows. Let the three tethers be indexed by a = 1,2,3. Each tether has a 
natural length L t a beyond which, when stretched, a restorative elastic force results. If the 
length of the tether is less than this, the tether will slacken and not exert any force in the 
opposite direction. Let the tethers have elastic constants given by k a . If the distance between 
the endpoints of the tether is given by L a , the magnitude of the elastic force due to tether a is 
given by: 

|^a| — k a AL a H a 
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where 

AL a — L a Lta 
and 

fl if AL a > 0 
a (0 if AL a < 0 

To determine the directions in which each force acts, we consider the case in which the tethers 
are between only two coils. The tethers are attached to individual points on the coils. For tether 

a, let the vector describing the point of attachment on coil 1 be given by P a , while the point of 

-+( 2 ) 

attachment on coil 2 is given by P0 . Thus, it is clear that: 




while the unit vector P e i as t,a parallel to the elastic force due to tether a is given by: 


5(D 


5(2) 


elast.a |-(i)_-(2)| 


With these definitions, we can write 


kaAL a H a P e ia S t,a 
ka.AL a H a P e ia S t,a 

where F a is the force due to tether a exerted on coil 1 and F a is the force due to tether a 
exerted on coil 2. The total force on each coil is then simply 

3 

£i(i) _ V 

r Elast / r a 
a= 1 

0 ( 2 ) 0 ( 1 ) 

r Elast r Elast 
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Now, the tethers will also effect torques upon the coils. If we wish to determine the torques 
about centers of mass of each coil due to tether a, we simply take the moment of the elastic 
force with respect to the coil center of mass: 


f (1) 

O, a, Elast 
-( 2 ) 

L 0,a,Elast 


— k a AL a H a 

— k a AL a H a 


(4 (1) - n) 

(p« <2) - r 2 ) 


* Pelast.a 
^ Pelast,a 


Thus, the total torque exerted by the tethers on each coil is given by 


= y f CD 
Zj °’ a ‘ 

a= 1 

= f f (2 > 
Z °’ a ‘ 


Elast 


r (1) 

O, Elast 

s(2) 

L 0, Elast / L 0, a, Elast 

a = 1 


One concern with this deployment strategy is that oscillations of the coils will occur around the 
equilibrium length of the tether. Potential solutions to avoid or damp out oscillations include a 
current controller that reduces and even switches the direction of current in the coils such that 
the coils experience an attractive force (i.e. a force in the opposite direction of their 
deployment motion) and are thus slowed to a stop as they reach the equilibrium length L eq . In 
this report, we do not explore a controlled method of damping in detail. 

Because any realistic tether will have intrinsic to it some damping due to friction between its 
constitutive materials. To account for this most generally, we would need to solve, 
simultaneously with those of the coils, the equations of motion of the tethers themselves. 
However, this is computationally intensive, requiring the additional solution of the three 
coupled partial differential equations (complete with tether-coil boundary conditions) of the 
elastic behavior of the tethers. 

For simplicity, therefore, we use instead a model in which the damping forces exerted on the 
coils are a function of the velocities of the tether attachment points only, as in a classical mass- 
damper system. Thus, we consider the relative velocity V a of the two coil attachment points for 
tether a with respect to the inertial coordinate system: 
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The damping force is supposed proportional to this, the proportionality constant taken to be 
b a . It can be shown, then, that: 


?(D 

a, damp 
?(2) 


\V a \H a 

Si g n (?a) a 

^ elast.a 

\Va\H a 

~ si g n (/a) ft 

* elast.a 


elast.a 


elast,a 


The total damping force acting on the coils is given by: 


= y f ( i) 

/ a, do 


' (1) - , , 

damp / a, damp 

a= 1 

F ( 2 ) = _p(l) 

r damp r damp 


The torques about the coil centers of mass due to damping from tether a are 

sign(? a ) 


f (1) 

B, a, Damp 
-( 2 ) 

L B, a, Damp 


| Va | 


ba | Kz | H a 


V3 

sign(F a ) 


V3 


' ^elast.a 


Pelast.a 


{Pa ] ~ fi) X P eiast , a ] 

(^a ^ — ^2) x A>iast,a] 


where the 6 denotes that the torques are taken with respect to the center of mass of the coil 
body. The total damping torques given by: 


- y f (i) 

Zj B ’° ’’ 

a= 1 

= f ? (2 > 

Zj B ’ a ’ 


Damp 


r (1) 

B.Damp 
s(2) 

L B.Damp / L B, a, Damp 
a= 1 
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Other options besides simple cable tethers attached at each end to a coil exist for the Separate 
deployment. Thin, rigid members with joints that lock open as the coils move apart in order to 
guarantee a fixed distance between the coils and reduce oscillations during deployment are 
also a possibility, though such constraints would be more massive and complex than simple 
tethers and thus less competitive in a design trade where mass is a driving parameter. The 
Separate coil configuration as modeled here is both a deployment and a support device; with 
rigid, locking members, it becomes purely a deployment device and thus loses the novelty of 
being able to perform two functions simultaneously. Another option, one that allows for 
Reconfiguration over the lifetime of the spacecraft, is tethers that can be reeled out from a 
spool or back in to change the unstretched length L t . 
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3.2.2 Hinges 

The other constrained structural configuration we shall consider is that in which two hinged 
panels are deployed by means of coil-generated forces. For the Unfold configuration, coils are 
embedded or otherwise attached to rigid panels aligned on one edge and connected with a 
hinge, with motion limited to rotation around that hinge. This configuration is useful for the 
deployment of solar panels or electronically steered radar arrays or other segmented flat 
structures. When current is run through the coils, the force that they exert on each other is 
converted to torque around the hinge, causing them to rotate towards one another for two 
currents in the same direction and away from each other for two currents in the opposite 
direction, as in Figure 13. 


In order for a hinged structure to deploy and maintain a predetermined equilibrium position, 
there must be a restorative force that opposes the coil driven separation. In practice, this will 
be accomplished by incorporating a torque spring or springs into the hinge about which the 
panels rotate. 

Flaving a spring also supports Reconfiguration or Deformation of a structure operationally, since 
shutting off current through the coils allows the spring to pull the coils back together. 
Reconfiguration or Deformation is also possible by reversing the current in one coil in order to 
cause attraction between the coils. 

Let us consider the torques about the hinge more in detail. Let the angles of the two panels 
around the hinge axis be given by 6 X and 0 2 . In the simplest case, if a torque spring is attached 
between the two panels, the spring torque is given by: 



Figure 13: Two coils repelling each other across a hinge in Unfold configuration 
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where k is the torque spring constant. Ad — 6 — 6 0 , 9 = \9 1 — 0 2 | is the angle between the 
two panels, and 9 0 is the natural angle of the spring. These angles are depicted in Figure 14, 
with the panels to which the coils are attached represented by lines from the bases of the coils 
to the hinge. 



Figure 14: Hinge at equilibrium angle of spring and at stretched angle 

Finally, in order to more effectively dampen any oscillations that occur about the equilibrium 
position, a rotational damping torque proportional to the angular velocity of the panels may be 
built into the hinge. 


In considering more fully the dynamics of the coils about a hinge, we first make some 
assumptions: 

1) The system comprises two panels, rotating about a hinge and in planes parallel to one 
another, with coils attached to each panel 

2) The hinge is fixed with respect to some inertial coordinate system 

3) The hinge is parallel to a unit vector which we denote j 

4) The coil is rigidly connected to the panel 

5) The coils and the panels are all rigid members (a brief discussion of what happens when 
the panels are not completely rigid may be found below) 
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These assumptions allow us to directly formulate our equations of motion. Since the hinge is 
assumed fixed with respect to the inertial coordinate system, we may use the parallel axis 
theorem in determining the moment of inertia of the system (since there is only one degree of 
freedom per coil in this system, we do not need to calculate the entire moment of inertia 
matrix). Let the distance from the hinge axis to the center of each coil as mounted on the panel 
be given by L m i , and let the radius of each coil be given by R t , i = 1,2. A coil rotating about an 
axis that passes through its diameter has a moment of inertia given by m t R 2 /2. The moment 
arm is given by L m i , so that, by the parallel axis theorem, the moment of inertia of the coil 
about the hinge axis is given by: 


r(0 

‘coil, Hinge 


m;R f , 

— + m iLm.i 


To get the moment of inertia of a single coil-panel construct, we must add to this the moment 
of inertia of the panel. Suppose the panel is rigid, rectangular, very thin compared to its length, 
and of mass 7n pane ^. If it has a length of h ir 

.(i) _ mpanel.ih-i 

1 P, Hinge ~ o 


so that the moment of inertia I hi f ge 


about the hinge for a single coil-panel construct is: 


Ji) m i^i . 2 m panel,i\ 

Hinge = + m i L rn,i + 3 


Again, the dynamics of the hinge deployment as specified allow for one degree of freedom of 
motion for each of the two panels, parameterized as the angle of each panel with respect to 
some axis. Then, Newton's 2 nd law can be written in angular form: 


/ 


d 2 0 1 ^ 

Hinge = 


( 1 ) 


/ 


( 2 ) d 2 9 2 

Hinge df2 


J = 


I 

I 


f (l) 

Hinge , 1 


-(2) 

L Hinge 


It remains to specify the torques acting on each panel. We will first consider the restorative 
torque applied to each panel counter to the moments due to electromagnetic force. These will 
be applied by torque springs about the hinge axis, and they may be implemented in one of two 
ways. 
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As mentioned above, the first way in which the torque spring may be implemented is by fixing it 
to each of the panels, so that the spring torque is a function of the difference in angle between 
the two panels, 9 — \9 2 — 9 ± \. Supposing the spring has a spring constant given by k and an 
unstretched angle between its end connections of 9 0 , the torque applied to each panel due to 
the spring will be given by: 


One consequence of this method of implementing a torque spring is that if there is any kind of 
jerking motion of either coil, then the coils will both acquire a net rotation, executing a motion 
in which the line bisecting the two panels has a nonzero angular rate about the hinge axis (as in 
Figure 14). This is undesirable when deploying the panels to precise orientations with respect to 
inertial space is required. Nevertheless, it is favorable from the point of view of minimizing the 
number of possible variables in the system and in that it directly ties the motions of the two 
panels to one another, if the deployment action requires so. 

As an alternative, we may consider the case in which each panel is attached to its own spring, 
each spring in turn connected to a structure that may be considered fixed. In this case, we have 
a spring constant for each panel, /q and k 2 , and an unstretched angle for each, 9 01 and 9 0 2 , so 
that the spring torques are given by 


The spring torques are now applied independently, but at the cost of having to determine 
additional parameters in order to implement a given equilibrium configuration. In this report, 
we will focus solely on implementing the first, connected, spring torques. 

As mentioned, we may also incorporate damping torques into the motion of the panels. How 
we do so, however, depends on the nature of the hinge itself. If the hinge is made up of two 
separate elements, each attached to a panel and rotating about an axle not connected to 
either, then the damping force will be proportional to the angular velocity of each panel 
independently about the axis: 





(Separate hinge elements) 
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where /? is the damping constant for the torque damping. If, however, the hinge is such that 
one element rotates about the other, then the damping torque will be a function of the 
difference in angular rates of each panel: 


Aamp= KK 
= -K(>i 


<y/i 

<y/j 


(Connected hinge elements) 


In this report we consider the hinge assembly to be made of separate hinge elements as in the 
first formulation. 

It remains to determine what the torques about the hinge axis are due to electromagnetic 
interactions between the coils. If we designate the meeting point of the panels and the hinge as 
the y-axis through the origin of the inertial coordinate system, then we have 


_( 1 ) _ Phff, 1 l eff , 2 

^ Hinge 4 n 


f t X dl t X 


Loop 1 


Loop 2 


dl 2 x r 12 
|r 12 l 3 


->( 2 ) _ dhff, 1 2 

^ Hinge 4 ^ 


r 2 


Loop 2 


x dl-, x 


Loop 1 


dli * ^21 

|r 21 l 3 


where r x is the distance from the origin to a point on coil 1, r 2 is the distance from the origin to 
a point on coil 2, and r 12 = r 2 — fy. The resulting torque will be in a direction parallel to the 
hinge axis. 
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3.3 Dynamic model implementation 

The equations of motion shown above were implemented and validated in MATLAB using a 
combination of built-in MATLAB functions, MathWorks' Simulink tool, and a number of custom- 
written scripts. 

Generally, dynamics involving the full Biot-Savart force law were computed in Simulink (which 
allows considerably for the incorporation of custom scripts), while those implementing the 
dipole force law were computed using the built-in ode45 solver. In the former case, the need 
to numerically compute the line integrals and resulting forces and torques about the coils 
necessitated the development of custom scripts, while in the latter case, analytical forms for 
the force and moment equations lent themselves to calculation using one of MATLAB's robust 
ordinary differential equation solvers. 


3.3.1 General solution algorithm 


The two-coil model was constructed for use in studying the actuation of circular, magnetic coils 
modeled as rigid, or unchanged from their circular state in the process of actuation. This model 
applies to the Unfold, and Separate deployment configurations and the Deform, Reconfigure, 
and Refocus operational configurations, all shown in Figure 15. Though this report focuses on 
two-coil systems, the model could be expanded to include three or more coils, which would 
enable study of the Inflate configuration. 



Figure 15: Rigid coil functional configurations 
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Below, we detail the steps taken in solution of the general equations of motion: 

• The model begins by defining the positions of cp evenly spaced points around each of 
two coils with respect to a global coordinate origin and assigns values for current (/ 12 ), 
wire cross-sectional area, wire density per unit length p, radii of the coils (f? 12 ), and the 
moments of inertia of each ring ([/] lj2 ). 

• At each timestep, the force and torque on each coil as a result of the magnetic field of 
the other coil must be calculated. This is done using a double for-loop that calculates 
the force and torque on each point around Coil 1 from every point on Coil 2 (and vice 
versa), summing the contributions on one point from Coil 2 and then summing all the 
resultant forces and torques around Coil 1 for a final force and torque value on the coil 
at that timestep. The resultant forces and torques are multiplied after the loops end by 
the currents and appropriate constants to complete the Biot-Savart equation and 
calculation of Laplace force. 

• Accelerations of the centers of mass are calculated by dividing the resultant 

— > 

electromagnetic force (F 12 ),) by the masses of the coils. 

• The positions of the coils are updated using the position and the velocity from the 
previous timestep. 

• To determine the rotation on the coil about its center of mass, the torques (one for each 
coil) are projected onto the local axes { X B } of the coils. The angular velocities, 
quaternions and direction cosines matrices are calculated from the torques, and the 
orientation of the coils in space are updated accordingly. 

• Plot new location and repeat in next time step. 

3.3.2 Note on solution of stiff equations 

-4 

F e iast is dependent on A L and the time is discretized by timesteps of varying duration in the 
tether model discussed above. We see from the formula for F eiast that the larger k is, the 

—4 —4 

larger the restorative F eiast force is for a certain A L. F eiast and the electromagnetic force 

— 4 — 4 

F magnetic act simultaneously over the course of a timestep At; as F magnetic stretches the 

— 4 — 4 

tether, F eiast is opposing this motion. The formula for F coii when L > L t is very stiff when using 

—4 

most materials. The F magnetic acting to accelerate the coils apart over a timestep At creates a 

—4 —4 

A L that causes a F eiast that is much larger than the F magnetic for all but the smallest ALs (and 
thus the smallest Ats are required for these two opposing forces to be of similar magnitudes 
and allow the model to converge to the approximation of continuous motion). A Simulink 
model utilizing the variable-timestep ode45 is used to reach the small Ats when L> L t and use 
larger Ats when L < L t for shorter overall run times. 
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3.4 Validation of numerical models 

Validation of the numerical methods employed in solving the dynamical equations of the coils 
are carried out by comparing results in which the force law is governed by the Biot-Savart law 
with those in the dipole approximation. 

Validation of the numerical approximation can be done via an approximation of the coils as 
magnetic dipoles at far-field distances. The far-field as discussed in this report is defined as the 
case in which the diameter D of each coil under consideration is much less than the distance 
between every pair of coils in the network: 


D i,j«Vij\. Uj = 1,2, 


A more thorough discussion of what constitutes this limit follows below. While validation 
against an exact analytic solution is not possible in the near field, we demonstrate that both the 
behavior of the numerical model in the near field is consistent with its behavior in the far-field, 
and that the behavior of the numerical model in the far-field matches that of the dipole. Thus, 
the validation provides confidence in the accuracy of the model. 

3.4.1 Analytic dipole model 

Now, the forces exerted by two magnetic dipoles on each other are given by Schweighart [20]: 


3p 

mi ■ rhj 

An 

r 5 


r + 




m 


— m j + 


■ f (mi ■ r)(m, ■ r) 

— T? 


where r is the vector distance between the centers of the coils, p is the magnetic constant, and 
in! and m 2 are the magnetic moments of each coil, defined by: 

ffii = UifAiUi 

— Ieff,i A-ifli 

n t being the number of turns in coil /, I t being the current, A t the area enclosed by the loop, and 
n i being the unit vector along the axis of the coil, oriented in the direction dictated by the right- 
hand rule following the current around the loop. 

For two identical, coaxially-aligned coils with planes parallel to one another, the dipole force 
equation reduces to: 
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-> A 1 A 2 I e ff ' 1 hff,2 „ 

f = s — fcr — ri2 

with f 12 the unit vector in the direction of r 12 towards the other coil. The above formula is 
multiplied by -1 if the currents are in opposite directions. 

3.4.2 Sources of error in numerical approximation 

There are two sources of error that contribute to the difference in force calculations between 
the numerical model and the far-field dipole approximation: 1) the discretization of the coils 
into straight line segments rather than a continuous curve, and 2) the distance between the 
coils. The errors resulting from these two sources decrease to zero as the number of line 
segments in the discretization or the distance between coils increases to infinity. 

Discretization: When using the Biot-Savart force law, the number of calculations to determine 
the resultant force and torque on each coil due to the other's magnetic field is 0 (<fi 2 ), with <p 
being the number of differential elements around the coil, making fine discretization of curved 
coils (though more accurate to the curve of the coil) computationally expensive when 
computing the electromagnetic force. 

Distance: The numerical force falls off as 1/r 2 , and the dipole force falls off as 1/r 4 , so we 
expect for very small (r < 1) values that the dipole force dominates, but then the dipole force 
is much greater than the numerical force right after r-1 and then the difference between the 
two calculated forces goes to zero as r -> oo. A quick visualization of this trend is shown in 
Figure 16. 



Figure 16: Proportional difference between numerical and dipole approximations 
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There is not a specific ratio of ; to r £j - at which the far-field approximation becomes valid; 
rather, one can determine what constitutes an acceptable percent error between the two. 
Figure 17 plots the initial force between a lm diameter coil and a 0.98m diameter coil facing 
each other against the distance between them in both the numerical model and the dipole 
approximation, showing the convergence of the initial force values as the distance between the 
coils increases. The coils are slightly different sizes in order to avoid overlap of the coils in 
space for very small (r £; - « 1) distances. 


Initial force between two coils for different starting distances 



Figure 17: Comparison of numerical and dipole initial forces over a range of starting distances 

Table 4 lists the error of the dipole force calculation compared to the numerical force 
calculation for distances that are various multiples of the coil radius. 

Table 4: Error between numerical and dipole force calculations 


Distance in coil radii 

Error (%) 

1 

87.4 

5.021 

16.5 

10.05 

4.7 

15.07 

2.2 

15.98 

2.0 

16.98 

1.8 

20 

1.3 
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Assuming an acceptable error in force of 2%, we find from Figure 17 and Table 4 that the 
distance at which the far-field approximation becomes valid is ~16 times the coil radius, or 8 D. 
This activity validates that the six degree-of-freedom approach presented in this chapter 
converges to the dipole approximation as expected at distances far from each other. There is 
no equivalent validation of the approach for the near-field, so this far-field convergence must 
suffice. We emphasize that the purpose of validation is to verify the accuracy and performance 
of our numerical methods in calculating the Biot-Savart force (which is a well-known 
consequence of the Lorentz force law and Gauss's law and is valid at any distance between 
coils) and to make sure that it agrees with another valid calculation of electromagnetic force 
that does not require a numerical approximation to calculate, the dipole force. The fact that 
the dipole case only converges to the true Biot-Savart forces in the far-field does not affect the 
Biot-Savart formulation in the near field: if the models agree in the far-field, then the Biot- 
Savart formulation is correctly and closely approximating behavior in the near-field. 

3.4.3 Validation models 

Validation of the numerical code for coil dynamics consisted of comparing the simulation 
results for full dynamics (i.e. those governed by the Biot-Savart force law) with those in the 
dipole approximation in the appropriate distance limit. Ten different deployment 
configurations were used for validation purposes: 1) Axially-aligned coils with elastic tethers 
between them, 2) hinged coil-panel systems with torque springs attached. For axially-aligned 
coils, rigid body rotation is not expected to occur. Therefore, validation for the Separate 
configuration focused solely on the rectilinear motion of the coils. The equations of motion 
implemented in the validation models will be discussed below. 


3. 4.3.1 Separate configuration 

For two axially aligned coils whose planes are parallel to one another, motion will exist (for 
either force model, dipole or numerical) in one direction. Thus, we may directly compare the 
axial separation of the two coils as a function of time; this will later serve to give an idea of 
what characteristic time scales are required deployment for various tether and coil sizes and, 
ultimately, for various load sizes. 

The equations of motion in the Laplace force case are: 


m l ^£2 ^ [ k a AL a H a P e iast,a + + 


P-hff, 1 hff, 2 
47T 


a= 1 

m 2 = ^ \kaAL a HaPelast,a + 4,damp] 
a= 1 


dli X 


Loop 1 Loop2 


dl 2 x r 12 

|r 12 l 3 


P-hff.l leff , 2 
47T 


din x 


Loop 2 Loop 1 


dli * r 21 

|r 2 il 3 


while in the dipole case, they are: 
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m 


1 dt 2 ~ [ kaAL a H a P e i ast a + F adam ^ + ^ 
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3. 4.3. 2 Unfold configuration 

In the hinge configuration, the equations of motion follow from the discussion above; for the 
Biot-Savart force law, we have: 


'HLsefjp-j = k(I »2 - 9,1 - e„); - tsej+j 


pJeff, 1 Jeff, 2 


An 


i\ x dli x 


dl 2 x r 12 


Loop 1 


Loop 2 


l?12l ; 


= -k(I»2 - 9,1 - «„)/ - 0J+1 


pJeff, 1 Jeff, 2 


4n 


f 2 x dl 2 x 


dli * r 2i 


Loop 2 


Loop 1 


|r 2 il : 


where I^i nge and J^ nge are defined as above. For symmetry and simplicity, the single torque 
spring model has been used. 


For the dipole case, we consider only the situation in which each hinged panel is of the same 
length. Consider the expression for the force acting on the centers of each coil, again letting 
A6 = \Q X — 0 2 1 be the angle between the hinge panels, we have by the geometry: 



3 ^ Jeff, 1 Jeff , 2 ^ 1^2 
An |r 12 | 4 


[1 + cos 2 (Z!0/2)], 


the force acting along a line that passes through the centers of the coils. If we calculate the 
torque due to this force about the hinge, we find (recalling that the moment arm is of length 

LmY- 


f (l) = 

Hinge 


Jeff, 1 
47T 


Jeff, 2 AiA 2 L r 


'12 


j( 2) 

1 Hinge 


[1 + cos 2 (Zl0/2)] cos(40/2)/ 
-f (1) 

Hinge 


Noting that r 12 = 2L m sin(40 / 2 ), we have, finally: 
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and 


f (1) = 

Hinge 54 ^ 

f (2) = -f (l) 

L Hinge L Hinge 


1 3/1 I e ff , i I e ff, 2 A\A 2 

1 + cos 2 (40/2) 

' 647T L m 4 

sin 4 (40/2) 

3 P Ieff, 1 Ieff , 2 A\A 2 

1 + cos 2 (40/2) 


sin 3 (zl0/2) tan(40/2) 
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Thus, the corresponding dipole model equations of motion are: 
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3p I eff, 1 0//,2 ^ 1^2 
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3p Ieff, 1 ^eff, 2 A\A 2 
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1 + cos 2 (40/2) 


sin 3 (40/2) tan(40/2) 
1 + cos 2 (40/2) 


sin 3 (40/2) tan(40/2) 


J 
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3.4.4 Validation results 

The foregoing equations of motion for the Separate and Unfold configurations were solved for 
several variations of the physical parameters, and the dipole and numerical results compared. 

In the Separate case, the coils are both of radius 0.5 m. They begin, axially aligned, at a distance 
of 5.0 m, with the tethers having an unstretched length L t — 6 m. The coil wires, assumed to 
have the density of copper, have a cross-sectional area of 10 -4 m x 4 x 10~ 3 m = 4x 
10~ 7 m 2 , have 90 turns, and current of 100 A running through them, which are typical values 
for industry standard HTS wire. [6] The baseline damping coefficient for a single tether is 
b 0 — 2/3 kg/s and the baseline single tether spring constant is k 0 — 0.1 N/m. The mass of a 
payload attached to each coil is 10 kg, which could include the weight of a cooling system for 
the coil or other apparatus. 

In the Unfold configuration, the coils have the same radius, cross sectional area, number of 
turns, and current as in the Separate configuration. The baseline damping coefficient 
/? 0 = 60 N • m • s/rad and the baseline torque spring constant is k 0 — 1 N • m/rad. The 
moment arm of the coil L m — 10 m, while the panel has a length of h — 20 m and a mass of 
m panei = 10 kg. The spring has an unstretched angle of 0 O = 45° = n/4 rad, and the coils 
have an initial separation of 9 init — 45° = 7t/ 4 rad. In both cases, the coils are discretized by 
0 = 17 points; these baseline parameters are summarized in Table 5. 


48 


Massachusetts Institute of Technology 


Space Systems Laboratory 


Electromagnetic deployment of space structures 


NASA NIAC Phase I Final Report 


Table 5: Baseline values of physical parameters for simulation 


Variables 

Symbol 

Value 

General 

Diameter 

D 

1 m (=2 R) 

Current through wire 

1 

100 A 

# of turns 

n 

90 

Effective current 

! eff 

I xn — 9000 A 

Density per unit length of wire 

P 

10~ 4 m x 4x 10 _3 m x 894 kg/m 3 
= 3.58 kg/m 

Discretization of coil 

<P 

17 points 

Separate Configuration 

Initial distance between coil centers 

d-stcurt 

5 m 

Unstretched length of tethers 

Lt 

6 m 

Tether damping coefficient 

b 0 

2/3 kg/s 

Elasticity of tethers 

k 0 

0.1 N/m 

Mass, payload 

m payload 

10 kg 

Mass, coil 

m 

pnD + m payload = 10.20 kg 

Unfold Configuration 

Initial angle between coils 

dinit 

45° or 7 t/ 4 rad 

Unstretched angle of spring 

do 

45° or 7 t/ 4 rad 

Torque damping coefficient 

Po 

60 N • m • s/rad 

Torsional spring constant 

K 0 

1 N • m/rad 

Mass of panel attached to coil 

m-panel 

10 kg 

Length of panel 

h 

20 m 

Hinge moment arm length 

km 

10 m 
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3. 4.4.1 Tether validation 

Coil separation as a function of time, numerical vs. dipole, for various k and baseline b 



Figure 18: Coil separation as a function of time for Separate configuration. Spring constant is varied, while 

damping coefficient is kept at the baseline value. 

In Figure 18, we see the effect of varying the spring constant k across all three tethers on the 
separation of the coils in the Separate configuration for both the numerical, Biot-Savart force 
law (solid line) and dipole approximation (dotted line). The differences between the numerical 
model and dipole calculations are consistent with the small errors predicted above in Section 
3.4.2. For all cases, the agreement between the numerical force calculation and the dipole 
approximation responses is good, with the discrepancy primarily due to the fact that the dipole 
approximation tends to overestimate forces at shorter separations. Since the coils start ten coil 
radii apart, agreement is within a few percent. 

The black lines show the response of the system as a function of time for the baseline case. 
Initially, the coils start 5 m apart, and they are accelerated away from one another by 
electromagnetic repulsion. After a separation of 6 m, which is the natural length of the tethers, 
elastic forces come into effect, and we see the coils begin to get pulled back while still being 
repelled from one another by the electromagnetic force. Ultimately, the motion is damped and 
reaches steady state after about 90 s. 
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The blue and red lines represent the response of the system for k = .1 x k 0 and k — 10 x 
k 0 , respectively. For the lower k system, we see the diminished effect of the elastic forces 
manifest in both the lower oscillation amplitude during pullback, and in the larger steady state 
separation acquired by the system; steady state is reached at a later time than in the baseline 
case. For the larger k, we see increased oscillation amplitudes as expected, while again, steady 
state occurs at a later time than the baseline situation. We can therefore surmise that, for the 
given b — b 0 , the baseline k represents something closer to a critically damped system than for 
the other two ks considered. For the higher k case, the steady state separation is smaller than 
the baseline case because the tethers are able to resist electromagnetic repulsive forces more 
effectively. 


Coil separation as a function of time, numerical vs. dipole, for various b and baseline k 



Figure 19: Coil separation as a function of time for Separate configuration. Damping coefficient is varied, while 

spring constant is kept at the baseline value. 

In Figure 19, we see the response of the system with the spring constant for the tethers fixed at 
k — k 0 , while b is varied. The black line represents the response of the system for baseline k 
and b values; this is the same response as that shown for the baseline case in Figure 18. Since k 
is not changed, the equilibrium value about which oscillations can occur is the same for all b 
values. 
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We see, for diminished b (b — . 1 x b 0 ; blue lines), that oscillations due to elastic forces (once 
again, coming into play once the separation is greater than 6 m) are less effectively damped 
than in the baseline case. For the larger, b — 10 x b 0 case (red lines), we see that the system 
is close to being critically damped, if slightly overdamped. Nevertheless, the high b case 
reaches steady state slightly sooner than the baseline case. Again, qualitative agreement 
between the numerical and dipole cases is good. 
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3. 4.4.2 Hinge validation 


Angular separation as a function of time, numerical vs. dipole, for various k and baseline p 
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Figure 20: Angular separation between hinge elements for the Unfold configuration. Hinge spring constant is 
varied, while torque damping coefficient is kept at the baseline value. 

In Figure 20, we see the angular separation response for the Unfold configuration as a function 
of time. Again, the black lines are the response of the system for the baseline case, for both 
numerical (solid line) and dipole (dotted line) force calculations. We see that the coils execute 
about two cycles of oscillation before equilibrium is reached at an angle of about 48°. For a 
lower torque spring constant k — .1 x k 0 (blue lines), the system is comparatively 
overdamped for the given torque damping coefficient /? = /? 0 . Subsequently, oscillations are 
nonexistent over the time period considered and steady state about a larger equilibrium 
separation is not acquired. On the other hand, for a larger torque spring constant k — 10 x 
fc 0 (red lines), we see a smaller equilibrium angle, due to the increased resistive capacity of the 
spring, and an increased number of oscillations before steady state is reached. 
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Angular separation as a function of time, numerical vs. dipole, for various p and baseline k 



Figure 21: Angular separation between hinge elements for the Unfold configuration. Torque damping coefficient 
is varied, while hinge spring constant is held at the baseline value. 

In Figure 21, we have the Unfold angular separation response for various torque damping 
coefficients /?, with the torque spring constant fixed at k — k 0 . Since k is fixed, the steady state 
angular separation of the system does not change, so that motions are either about this steady 
state value (equal again to about 48°) or asymptotic to it in an overdamped case. The baseline 
case (black) is the same as in Figure 20, while the response for a decreased damping coefficient 
/? = .1 x /? 0 is shown in blue. We see an increase in oscillation cycles for a given time period 
over the baseline case due to the system being relatively underdamped, while by comparison, 
the /? = 10 x /? 0 case is overdamped (red lines), not reaching the steady state angular 
separation and instead asymptotic to it. 
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3. 4.4.3 Effect of discretization 

Coil separation as a function of time, numerical vs. dipole, for various coil discretizations tj) 



Figure 22: Effect of discretization on numerical calculations. Baseline tether case is calculated for Biot- 
Savart/numerical and dipole force laws, with the former approximated by increasingly finer discretizations 0. 

The inset shows the same, near the oscillation peak occurring around 20 s. 

In the previous sections, we investigated the response of the Separate and Unfold 
configurations as a function of time for both the dipole and Biot-Savart force laws. The 
numerical calculation of the latter depended upon the discretization 0 of the coil, with 0 the 
number of discrete points approximating the circular shape of the coil. Again, the baseline 0 
used in the previous calculations was 17. Nevertheless, it is important to understand the effect 
of varying 0 from the standpoint of numerical accuracy, while for computational efficiency, it is 
worthwhile to keep 0 as small as possible without sacrificing that accuracy. 

In Figure 22, we show the effect of varying 0 on the response of the Separate configuration, 
with baseline parameters used for k and b. For reference, the dipole calculation is shown as the 
dotted black line. The smallest discretization considered, 0 = 9, is shown in solid black, while 
the responses for 0 = 17, 25, and 33 are shown in blue, red, and green, respectively. While the 
error is considerable for the 0 = 9 case, the 0 = 17, 25, and 33 cases are close to one another, 
with the 0 = 17 and 25 approaching the 0 = 33 line, which is naturally closest to the true 
response. Since for each integration time step, the number of calculations goes as 0 2 , we have 
chosen 0 = 17 for both its relative speed and accuracy. 


55 


Massachusetts Institute of Technology 


Space Systems Laboratory 


Electromagnetic deployment of space structures 


NASA NIAC Phase I Final Report 


3.4.5 Note on elasticity of panel 

We note that there exists an expression for the torque about the center of mass of a coil in the 
dipole limit. The expression is 


While there will be a tendency for the magnetic moments to want to rotate about their centers 
due to the above torques, it is important to note that, as long as the panels are taken to be 
rigid, these torques exert no influence on the motion of the coil about the hinge axis. The 
influence of these torques would only be manifest if the panels were elastic, in which case the 
torques around the center of the coils would impart a moment on the panel resulting in 
bending of the panel. Accounting for this would necessitate solving the dynamic bending 
equations for the panels and would, to an extent determined by the deployment geometry and 
the material properties of the panel, affect the distance between the coils. 

To gauge what this criterion would be for which bending is appreciable, we may conduct a 
simple order of magnitude analysis. From the theory of beam bending, we have that the 
bending moment of a beam is given by: 


where x is a coordinate directed along the axis of the unbent beam, E is the elastic modulus of 
the beam material, and J (not to be confused with the current) is the area moment of inertia of 
the beam. If we approximate the deflection y of the beam due to the dipole torque as 8 and 
take x ~ L m , then we have, using the dipole torque expression. 





E hff, 1 hff, 2 A\A 2 


or, rearranging. 



2n r 3 EJ 
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Now, the criterion for which this would have a negligible effect on the hinge dynamics of the 
coil-panel system (by altering the coil center-to-coil center distance) is that: 1 

8 

- « 1 
r 


or 


u L1 2 A- [ A 2 L 2 
— — « 1 
2n r 4 EJ 

Since E is typically on the order of 10 10 — 10 11 Pa for metals and p = 1.26 x 10 -6 T-m /A, we 
can rewrite the above in a more illuminating manner by considering a lower elastic modulus 
case: 


hhA\A 2 
Jr 4 


A 2 

« 2n x 10 16 — 
m / 


This criterion will generally be easily satisfied. Furthermore, since at short distances, the dipole 
approximation for forces and torques tends to overestimate that due to the actual force law 
(see below), the above expression may be taken as a general rule for ensuring that bending 
effects are negligible. 


3.4.6 Motion of a coil under its own force in Expand configuration 


Consider a flexible coil that is in some noncircular configuration (as in Figure 23). If current is 
run through this coil, a given differential coil element will experience a force due to the 
electromagnetic forces exerted on it by all other elements on the coil. It is therefore possible 
for a stowed coil to reconfigure itself, or Expand, into a different shape under the influence of 
its own electromagnetic interactions. We note that, for the rigid coil considered elsewhere in 
this report, such self-interactions have no bearing on the shape of the coil, so that in that case, 
it is assumed that the coil can effectively "absorb" the effect of self-interactions. 



Figure 23: Expand deployment configuration 


1 A more precise statement may 


be 



'r 12 


r 12 


S 


« 1; however, for simplicity, we use the 


given relation. 
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Suppose we have a flexible coil, the position of a point on the coil given by r. Let us 
parameterize the position of a point on the coil by the coil length s, so that, in general, 
r = r(s, t), where t is time. If the string is of constant length (or inextensible), the equations of 
motion for the string are given by: 

d 2 f d f dr\ 

P dF = d^\ a d^) + F 

dr dr 
ds ds 


The second equation is a restatement of the Pythagorean theorem as applied to a infinitesimal 
coil element, o is the magnitude of the local tension in the string, so that the vector 

dr 
a — 


is a force acting across the cross section of the coil. For an inelastic coil, both equations above 
must be solved simultaneously. If however the coil is elastic, following a linear elastic force law, 
the equations of motion reduce to a generalized wave equation: 


d 2 f 

P di 2 


d 

ds 



+ F 


where E is the elastic modulus of the coil material and A c is the cross sectional area of the coil. 
In general, both are a function of s, but, if the coil material is uniform along its length, we have: 


d 2 f d 2 f 
PTTIT = ea c — + F 


dt 2 


ds 2 


— ) 

The force F is calculated from the Biot-Savart law given before, but applied to the coil itself: 


-> pi 2 dl f Idl' x (r — r') 

F 4tt Ufl X J \r — r'| 3 

1 l coil 
f=£ rr 

where f' is the vector to any point on the coil not equal to r and dl' is the differential coil 
element at f '. An in depth discussion of the dynamics of self-coil deployment will be reserved 
for the Phase II report; nevertheless, we include here some preliminary results of calculations 
relevant to the self-force deployment. 
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Figure 24: Instantaneous self force on coil for two different coil shapes. Direction and magnitude of force 

indicated in arbitrary units. 


In Figure 24, we show the instantaneous force vectors acting along coils of different shapes. 
We note that inflection points along the curve are coincident with locations where the 
electromagnetic self-force vanishes. 

In Figure 25, we see a visualization of the time evolution of self-deployment of a coil in 
accordance with the equations of motion above for an inelastic coil. The coil starts from a 
stowed, folded position, and, over time, acquires a circular shape due to the mutual repulsion 
of coil elements. We note that, during self-deployment, the coil will generally oscillate for 
some period about an equilibrium circular configuration. Nevertheless, as any material is not 
perfectly elastic, these oscillatory motions will ultimately be damped out, the efficiency with 
which is occurs depending on the construction and material properties of the coil. At the 
present, we have not implemented elasticity in the coil material. 
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Figure 25: Visualization of self-deployment of a coil from a stowed configuration to a planar, circular shape. 


3.5 Conclusion 

In this chapter, we have detailed the derivation of the equations of motion for current-carrying 
coils under the influence of electromagnetic forces and constraining elements. We have 
developed and validated a set of numerical simulation tools to evaluate the performance of 
electromagnetic system dynamics. By using the numerical models for rigid coils with 
appropriate mechanical constraints, we can now begin to set up the systems engineering trades 
for the implementation of electromagnetic deployables in candidate space system designs and 
investigate the ranges of values over which an electromagnetic structure is feasible. 
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4 Chapter 4 - Feasibility 

In this chapter, we will study the feasibility of electromagnetic structures in space systems by using 
the models described in Chapter 3 to generate data over a range of design variables for the 
Separate and Unfold functional configurations, including current, coil diameter, number of turns, 
and others. These data are intended to provide a means of trading electromagnetic designs with 
one another for different application requirements. We will also introduce the prevailing trades 
between variables and between electromagnetic structures and alternative structural technologies 
such as inflatables and tensegrity structures. 

4.1 Electromagnetic structure design variables 

Electromagnetic structures have a variety of independent design variables as well as dependent 
design variables that are a function of the other variables. Table 6 displays several of the most 
important variables in the design vector for electromagnetic structures as well as their 
interdependencies, which are listed in the third column. 


Table 6: Design variables 


Variables 

Symbol 

Function of? 

Independent, General 

Current through wire 

/ 

-- 

Diameter 

D 

-- 

# of turns 

n 

-- 

Density per unit length of wire 

P 

-- 

Initial distance between coil centers 

d-start 

-- 

Initial angle between coils 

@init 

-- 

Angle between coils (overtime) 

e 

- 

Stowed structure length 

L s 

-- 

Damping coefficient 

b,P 

-- 

Discretization of coil 

<P 

-- 

Independent, Tethered 

Unstretched length of tethers 

L t 

-- 

Elasticity of tethers 

k 

-- 

Independent, Hinged 

Torsional spring constant 

K 

-- 

Mass of panel attached to coil 

Impanel 

-- 

Unstretched angle of spring 

8o 

-- 
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Dependent 

Effective current 

hff 

\, n 

Current element 

Id! 

1 eff> V/ D 

Mass (each coil) 

m 

D, n,p 

Mass (total) 

™ tot 

TH-panel 

Distances from elements on coil 1 to coil 2 

r a 

D, d, &, (p 

Force (one element on another) 

F 

Idl, r tJ 

Force (total on coil) 

Ftot 

F 

Distance between coil centers (over time) 

d 

d 0 , F, K, HI, ITlpanel 

Force (dipole) 

Fdip 

I ef f, D, n, d, & 

Acceleration 

a 

F, m 

Characteristic/Equilibrium length 

^eq 

F, d, L t and k (if applicable) 

Size ratio 

s 

Leq* L s 

Torque 

T 

F,D 

Hinge moment arm length 

L-m 

d/ @init 

Equilibrium angle of hinge 

S eq 

F, k, 

Deployment time 

T 

F dip or F, m, L t , b, k,G 0 , L m , d 

Moment of inertia 

VI 

^hinge 

L-mi ®init' D, Impanel 


There are several other design variables that come into play when more than two coils are being 
considered, the coils have variable current controllers, or when the coils are different sizes, but for 
now we are just examining two-coil, same-sized coil systems with constant currents throughout the 
deployment. 
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4.2 Reformulation of equations of motion for trade space analysis 

Before exploring the variable or trade space, it is important to note that the number of 
independent variables, when the equations of motion are written as is, is large. To examine 
variance across each of the dependent variables as a function of changing each independent 
variable individually would be an overwhelming amount of data and plots. Therefore, it is 
worth exploring methods by which the dimensionality of the variable space can be reduced. 
One method by which this can be accomplished is via non-dimensionalization of the equations 
of motion. This process, which is commonly used in fluid dynamics and heat and mass transfer, 
greatly reduces the variable space by creating dimensionless variables that are combinations of 
the original ones. The process will be demonstrated below for both the Separate and Unfold 
configurations. We assume that in both cases the relevant geometries and properties for each 
coil, panel, and tether are the same (i.e. m x — m 2 = m, k 1 = k 2 = k, etc.). 

First, we choose an important, characteristic geometrical variable against which we scale all 
variables in the problem that have dimensions of length; next, we do the same for the time 
scale. We choose R, the radius of each coil, as the characteristic length variable. We choose as 
the characteristic time variable T, which we set arbitrarily to 1 second. Now, we proceed to 
define dimensionless auxiliary variables based on the variable originally appearing in the 
equations of motion; these dimensionless auxiliary variables will be distinguished by tildes 
or asterisks "*" and are given in Table 7. 
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Table 7: Dimensionless variable definitions 


Dimensionless variable 

Variable definition 

R 

Ri=- = 1 
1 R 

Radius of coil / 

« _ n 

n ~ R 

Vector position center of coil / 

* _ n 
n ~ R 

Vector position of a point on coil / 

5 _nj 

r i J - R 

Vector distance between a point on coil /' to a point on coil j 

^ 1°= 
hi 

Vector differential length element at a point on coil / 

. * _ b t?a 

L't.a — ^ 

Length of tether a 

v / 

III 

Vector position of tether a attachment point on coil i 

t 

£ = T 

Time 

Tts? 
fc-> 1 Oh 
III 

itts? 

Relative velocity of tether a attachment points on coils 1 
and 2 

j * _ d s tart 

“■start — R 

Starting coil separation 

, * _ ^eq 

eq ~ R 

Equilibrium separation of coils 
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Substituting the above dimensionless variables in for the original variables yields the following 
general equations of motion: 

Tether: 


d 2 (Rr 1 ) 
1 d(Tt) 2 


: - ^ k a (RAL a )H a P elast . a + b a \v a \ H a 


sign (l4) 


V3 


* plant, nr 1 




47T 


(Rdf) 


Loop 1 Loop 2 \^ r 12 


d 2 (Rr 2 ) 
m d(Tt) 2 

Hinge: 


: + ^ k a (RAL a )H a P elast , a + b a (£) \v a \ H a 


sign 


(l) 


V3 


f (Rdf), f 

Loop 2 Loop 1 (^^21 1 


CD d 2 8 1 

Hin a e d(Tty 


; = «CI S2 - Sll -e„ V -pf^ t + ^a^ 


Loop 1 


Loop 2 


C ge 4Jk;J = -^C| 6» 2 -6,1- Go)] - P + 


d(Tt) 2 


d (T t) An 


(Rrf) x (^Rdl^) x 
j (R?z) x (i?d? 2 ) x j 


(Rdl^) x (Rr 12 ) 


d0 2 f , lileff.ileff .2 X r n Z^fn, 7 \„ X 

Loop 1 1^*21 1 


Loop 2 


Collecting all the constant terms then yields: 
Tether: 


d 2 f 1 
~dF 


d 2 f 2 

TO 2 




sign 




V3 ‘ e!ast ’“ 

sign ( 14 ) 


Neff.l hff, 2 T 2 

An mR 


Loop 1 


V3 


Neff. 1 4 //, 2 
AnmR 


dl 2 x f 12 
|4 [3 

Loop 2 Xl2 | 

dll X ?21 


Loop 2 Loop 1 I ' 2 1 1 


Hinge: 


d 2 0i kT 2 BT d6, pI e fflIeff2‘J' 2 R 

-J = + (1 0 2 - e, | - 0 O )/ - T~ ST J + ff 


dt 2 


d 2 6 2 

dt 2 


I Hinge 

KT 2 

l Hinge 


I Hinge dt 


AnI. 


Hinge 


r, x dL x 


dl 2 X T 12 


Loop 1 


BT dO? B^eff 1 leff 2 

(102 - 01 1 - 0o); - T~ TrJ + 


I Hinge dt 


l Hinge 


f 2 X d/ 2 x 


Loop 2 


Loop 2 r 12 j 

dl 1 X r 2 i 

l^lT 


Loop 1 
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Next, we will redefine the dimensionless coefficients from the above equations in order to 
simplify the expressions. We also prescribe "baseline" values to each of the constituent 
variables in order to find a baseline value for each of the dimensionless parameters to use in 
anchoring plots of variable trades. The dimensionless quantities presented below in Table 8 
now represent the key parameters for each functional configuration: 


Table 8: Dimensionless parameter definitions 


Dimensionless 

parameter 

Parameter 

definition 

Baseline values of 
included variables 

Dimensionless 
parameter values for 
trades [Baseline] 

khff, 1 hff, 2 F 2 

v = 

4tt mR 

Tether (Separate) 
force coefficient 

hff , i > hff, 2 — 9000 A 
m = 10.2 kg 
R — 0.5 m 

y 0 = 1.588 

k a T 2 

k* = 

m 

Tether (Separate) 
spring constant 

k a = 0.1 N/m 
m — 10.2 kg 

k 0 * =0.010 

KT 

m 

Tether (Separate) 
damping constant 

b a = 2/3 N/m 2 
m — 10.2 kg 

b 0 * = 0.065 

c _ B-hff.l hff, 2 T 2 R 
4 ^-Itlinge 

Hinge (Unfold) 
force coefficient 

hff.l < hff, 2 — 9000 A 
R — 0.5 m 

L Hinge 

— 1.343 x 10 3 kg-m 2 

m 

o 

o 

II 

o 

k\/\ 

kT 2 

K* = 

1 Hinge 

Hinge (Unfold) 
spring constant 

k — 1 N • m/rad 

1 Hinge 

— 1.343 x 10 3 kg-m 2 

Ko * =7.443 X 10~ 4 

BT 

P*=T^~ 

l Hinge 

Hinge (Unfold) 
damping constant 

/? = 60 N • m • s/rad 

^ Hinge 

— 1.343 x 10 3 kg-m 2 

Po = 0.045 


We are able to plot dimensionless dependent variables of interest against these dimensionless 
parameters in order to study the effects of certain system trades in each of the Separate and 
Unfold configurations as well as the behavior of unconstrained coils. Section 4.3 introduces and 
models several of the more interesting trades. 
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4.3 Trades 

In each of the three models/functional configurations introduced in this report, there are a 
number of trades of system parameters that can be investigated in order to see relative 
sensitivity between variables as well as establish what combinations of parameters result in a 
viable system design. This section will discuss these trades from a mission level perspective and 
then explore them in detail using the models we have constructed. 

Three of the most important parameters across most space missions are mass, required power, 
and system cost, reduction being desired in each of the three. As can be seen from the 
definition of the coefficient y, force is proportional to the product of mass and current (power 
being proportional to the current through each coil as well as the number of coils in the system, 
which is fixed at two in this analysis). Therefore, mass and current can be traded against each 
other for a fixed y, or y can be reduced to a small value, thus minimizing both. The power of 
nondimensionalizing the variables in this problem is the ability to group these important 
parameters together such that the interior parameters can be switched around and still retain 
the same dimensionless parameter value and thus system behavior. A system designer can take 
the plots in this section and determine the desired y or f for the desired behavior, as well as 
corresponding spring and damping constants, and then use the defined relationship to trade 
the interior parameters like mass and current. System cost is assumed to correlate positively 
with the number of turns in each coil and with the number of coils in the system, so lowering 
mass and current can be related to a reduction in system cost. Note that a complete cost 
analysis was not performed in the course of this Phase I work but is planned for future 
hardware development research. 

Many of the plots in this section plot y or f against deployment time T for various k* and b*. 
While deployment time is not usually a key performance parameter in the deployment of space 
structures, these plots are mainly for exploring reasonable ranges of values for the system 
parameters. Plotting against deployment time gives a nice visualization of how changing each 
variable affects the end result of the system, as well as lets us see easily when a certain design 
does not ever converge to a steady state. Also, some of the operational configurations that are 
possible with electromagnetic structures may be more time critical, such as a Deformation to 
avoid a collision or a Reconfiguration quickly in reaction to a sudden event of interest for 
observational spacecraft, making deployment time an important parameter in these cases. 
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4.3.1 Six Degree of Freedom ( 6 D 0 F) free 

The 6D0F free case does not correspond to any of the functional configurations in particular, 
but it represents the basis for all electromagnetic structures and is valuable for studying some 
of the most basic relationships between variables without the complication of boundary 
constraints, especially between force and resultant motion. Force is very important at a 
mission level because all motion and equilibrium configurations are directly related to the force 
coefficient. 

The force generated between two electromagnetic coils is proportional to 1/d 2 By 
nondime nsionalizing d into multiples of coil radii (d/R) and using a range of ys as specified in Table 
8, we can see easily the effect on the force of variation of the variables contributing to the force 
calculation. Figure 26 relates the force coefficient y to the time it takes for two unconstrained coils 
to repel to twenty times their original starting distance. 



Force coefficient y 


Figure 26: Time T to separate to 20 x d start * as a function of y plotted for values of d start * from 0. 5 to 5. 5. 
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Deployment time decreases with increased force, which follows from Newton's second law. 
Since force is comprised of the force coefficient y and a term dependent on distance between 
the coils, the distance term dominates the force calculation for small ys, especially for small 
d-start *■$> an d we see that the lines are divergent and that 7” is influenced strongly by d start * for 
small ys. High ys dominate the force calculation and so we see convergence and little 
dependence by T on d start with high ys. 
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4.3.2 Tethered constraints in Separate configuration 

The tethered case presented in this Phase I report examines the behavior of axially aligned coils 
during separation, or the configuration pictured in Figure 10 with three tethers. The behavior 
of coils with altered angles (which introduces torques on the coils), axial offsets or attraction 
instead of repulsion will be investigated further in our Phase II work. 

4.3. 2.1 Force F tot versus deployment time T 

The deployment time of two tethered coils depends on the electromagnetic force between 
them, and the damping coefficient b* and spring constant k* of the tethers. The values of k* 
used in the plots in this chapter are for each tether; effectively, the spring constants are three 
times this value since there are three tethers in the system. The electromagnetic force 
depends on two parameters in the nondimensionalized equations of motion: y and the distance 
|fi — f 2 1 between the coils. Deployment time in the tethered case is defined as the time the 
system takes to reach the steady state. For numerical calculations, we define the steady state 
to be the case in which both the relative velocity and acceleration of the two coils are arbitrarily 
close to zero: 


d(r i - f 2 ) 


d 2 (r i -f 2 ) 

dt 

— ci, 

dt 2 


For our calculations, we choose: 

e 1( e 2 = 2x 10" 6 

The following plots investigate how the deployment time T is a function of y, b * , k*, d start *, 
and L t *. These variables encompass most of those variables from Table 6 that affect the 
motion and steady state position of the system without our plots becoming redundant with 
variables that are closely related to one another. In each of the following cases, one variable is 
being varied, and the others are fixed at the baseline values given in Table 8. 
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4.3. 2. 1.1 Varying damping coefficient b* 

Damping in our tether model occurs only while the tether is being stretched past its 
unstretched, natural length L t *. Damping depends on the tether materials and its construction. 
Oscillations are not desirable in deployment of space structures, especially large ones, because 
they could possibly amplify or aggravate structural incongruities or defects and the effect of 
such defects on stability of the structure. We can use deployment time as a metric to see 
quickly which system designs result in the shortest period of oscillations before convergence to 
steady state. Note that this plot does not tell us anything about the magnitude of the 
oscillations, just the duration. Section 3.4.4 shows validation plots of the motion of coils over 
time, highlighting the amplitude and duration of oscillations. Figure 27 shows how the 
deployment time T varies as a function of y for variable b*, with k* and d start * fixed at their 
baseline values. The baseline value, b 0 *, is drawn in red. 

Deployment time 7as a function of force coefficient y for various b* and baseline k* 



Figure 27: Deployment time T as a function of y plotted for different values of b* from 10 5 b 0 * to 10 b 0 *. 

We can see that T is higher (underdamped) across all ys for lower values of b* than for higher 
values of b* until we reach the minimum deployment times at 10 25 b 0 * (teal), which represents 
a value close to critical damping for the system. 10' 5 x b 0 * (purple) damps out oscillations 
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slightly faster than teal for high values of y. Then, starting with 10 5 x b 0 * (purple) and 
continuing on through the progressively larger values of 10' 75 x b 0 * (gold) and 10 x b 0 * (grey), 
we see a trend back up towards longer, overdamped Ts. The general trend is towards longer Ts 
from smaller y s, since smaller electromagnetic forces result in smaller accelerations, and, 
hence, smaller damping forces and longer oscillatory motion. For the values of damping at and 
around the baseline (red, teal), the deployment time remains fairly consistent across the span 
of ys (after an initial startup lag due to the small forces) because larger electromagnetic forces 
are met with larger damping forces in response to their higher velocity and are damped out just 
as quickly as medium forces in this critical damping regime. 
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4.3. 2. 1.2 Varying spring constant k* 

Figure 28 shows how the deployment time T varies as a function of y for variable k* and with 
b * and d 0 * fixed at their baseline values. The baseline value k 0 * is drawn in red. For reference 
and comparison with Figure 27, the red baseline line is the same in each plot. There is a clearer 
order in Figure 28 than in Figure 27: small ys lead to longer Ts and larger ys lead to shorter Ts. 
This is in accord with the observation that smaller electromagnetic forces, for a given b*, give 
rise to smaller accelerations and hence smaller damping and longer periods of oscillation. 


Deployment time Tas a function of force 
coefficient yfor various k* and baseline b* 



Figure 28: Deployment time T as a function of y plotted for different values of k* from 10 s k 0 * to 10fc 0 *. 
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The baseline (red) and the two smaller-than-baseline k* values (blue and green) give the 
shortest Ts for all ys plotted. We see a larger spread of Ts between the k* values with small ys, 
with larger k*s leading to longer T. Nevertheless, all k* lines converge to nearly the same T for 
large ys. This narrowing follows from the increase in damping due to larger electromagnetic 
forces. Coils moving initially with the same small y (and thus the same velocity) are pulled back 
together much harder upon reaching the stretching point of tethers with high spring constants. 
Large spring forces increase the number of oscillations before settling to the steady state 
length, thus increasing T. With large ys, the large resultant velocity of the coils triggers a 
correspondingly large damping force. Larger k*s and their larger resultant spring force and 
snapback velocity are also acted upon by damping, so even though they would, sans damping, 
take much longer to settle to steady state than cases with small k* s, damping serves to balance 
that effect. 
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4.3. 2. 1.3 Varying initial distance d start * 

Figure 29 shows how the deployment time T varies as a function of y for variable d start *, with 
b * and k* fixed at their baseline values. The baseline value, d start * = 10 (coil radii), is drawn in 
gold. For reference and comparison with Figure 27 and Figure 28, the gold baseline line is the 
same as the red baseline line in the previous two plots. Usually when deploying from a stowed 
position, the starting distance between the coils is small (like the blue line), but studying how 
starting position affects the deployment behavior lets us consider multi-stage deployments 
where the coils start farther apart after an initial deployment activity. We use 10 coil radii for a 
baseline starting distance throughout this chapter arbitrarily, but Figure 29 conveys the effect 
of changing that value on the system dynamics. 


Deployment time las a function of force coefficient y 
for various d , * and baseline b*, k* 

start 
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The behavior in Figure 29 is similar to that in Figure 28; small ys lead to longer Ts, and larger y s 
lead to shorter Ts, but the various d start * values converge to the same value of T at y ~ 12 in 
Figure 28 while they have nowhere near converged by y ~ 100 in Figure 29. The range of Ts 
for y ~ 0.1 is wider in Figure 28, however. This observation means that for small y < 0.5, the 
system is more sensitive to changes in k* than in d start *, but for y > 12, the system is almost 
completely insensitive to changes in k*, while it remains sensitive to changes in d start *. When 
the coils start closer to the unstretched length of the coil (in this case, L eq * = 6 coil radii), they 
have less time to accelerate before their motion is resisted by the tether spring force. That is 
why d start * = 5.5 (black) deploys the fastest regardless of y's value; the coils do not gain a high 
momentum that must then be damped out to achieve equilibrium. The coils that start closer 
together (small d start *s) experience exponentially higher electromagnetic forces than those 
coils that start farther apart due to the force's proportionality to the inverse of the distance 
between the coils squared. Larger forces result in larger oscillations that take more time to 
damp to steady state, so small d start *s have longer Ts across the board. 
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4.3. 2. 1.4 Varying tether length L t * 

In Figure 30 and Figure 31, we see how, in the Separate configuration, the deployment time 
varies as a function of the force coefficient y for various dimensionless unstretched tether 
lengths L t *. The baseline value, L t * = 12 (coil radii), is drawn in green. The green line is the 
same baseline line that is present in the other y versus T plots. The length of the tether is 
directly related to the characteristic or equilibrium length of the structure, L eq *, which will be 
discussed further in the next section. Since d start * is fixed at 10 coil radii in each case, we can 
see how increasing L t * for a given starting distance affects the deployment behavior. Figure 30 
and Figure 31 are very similar to Figure 29 because the difference between the starting and 
unstretched lengths is what really matters, not the value of either in particular. The values of 
each of L t * and d start * affect the component of the force that includes distance, but y can be 
adjusted to balance that effect, leaving the relationship between L t * and d start * and T 
primarily dependent upon L t * - d start *. 

Deployment time Tas a function of force coefficient y for various L* and baseline b*, k* 



Figure 30: Deployment time T as a function of y plotted for different values of L t * from 5 to 10 coil radii. 
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In Figure 30 and Figure 31, we see how, in the Separate configuration, the deployment time 
varies as a function of the force coefficient y for various dimensionless unstretched tether 
lengths L t *. The spring and damping coefficients k* and b*, and the starting coil separation 
d-start* are the same for all cases. For low y we see how, expectedly, deployment time is longer 
for larger L t *. This longer time is because a greater distance must be traversed from a given 
d-start* for the tether's restorative force to take effect. This is most clearly demonstrated in 
Figure 31, where long tether lengths of 15L t * and 20 L t * have deployment times well over 1000 
s (The model was only run for 1000 seconds to generate Figure 31, so the y-intersect for the 
two longest tether lengths are not included but can be extrapolated from the other points). On 
the other hand, for larger ys, we see that the acceleration of the coils (and subsequent 
damping to the steady state) occurs so fast that all the deployment times for different L t * tend 
to converge to the same T. 



Force coefficient y 

Figure 31: Deployment time T as a function of y plotted for different values of L t * from 5 to 20 coil radii. 
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4.3. 2. 2 Force F tot versus characteristic/equilibrium length L eq * 

The equilibrium length L eq * of the deployed structure is a very important metric for spacecraft 
designers, as high-level performance requirements generally dictate the deployed size. L eq * is 
thus more directly representative of the viability of HTS structures as compared to other 
structural technologies than is deployment time. The major design variable besides the force 
coefficient that affects L eq * is k*. Damping (b*) and the initial position of the coils d start *only 
affect the deployment time, not the steady state length, so long as the tethers are assumed to 
stay below their yield stress and therefore not stretch inelastically. Figure 32, Figure 33, and 
Figure 34 show these effects very succinctly. In Figure 32, where k* is varied, we see that 
different k* s result in different values of L eq *, as do different L t * s in Figure 33. Flowever, in 
Figure 34, which gives the same single line regardless of how many b* or d start * values are 
being plotted, we see that neither b* nor d start * affects the L eq * for a given y. The takeaway 
from this trade is that only the force coefficient as well as the spring constant k* and natural, 
unstretched length L t * of the tethers have the ability to affect the final, steady state 
configuration of the tethered structure. This is consistent with the conceptual model, because 
these parameters define the force balance in the steady state condition, while other 
parameters only define the starting conditions or describe the intervening system dynamics 
during deployment. 


4.3. 2. 2.1 Varying k* 

Figure 32 shows how the equilibrium length L eq * varies as a function of y for variable k*, with 
b* and d 0 * fixed at their baseline values. The baseline value of k* is the red line. Figure 32, 
Figure 28, and Figure 31 highlight one of the more interesting system trades in a tethered 
system: in order to increase the equilibrium length of the system without adding coils, one can 
employ tethers with longer unstretched lengths or use tethers with shorter lengths that require 
stretching. Longer tether systems do not require as much current in each coil after initial 
deployment, since electromagnetic force required to maintain coils at the taut length of the 
tethers is significantly less than that required to actually stretch the tethers. 

However, as we see in Figure 31, shorter tethers take less time to deploy than longer tethers, 
and low k* values deploy the fastest, as we see in Figure 28, because damping occurs during 
periods of the tethers being stretched and low k* tethers stretch the most. We assume mass- 
less tethers in the models in this report, but shorter tethers also present mass savings, though 
on a scale that would only matter to relatively small spacecraft. Low k* tethers also provide the 
capability for a system to Reconfigure or for some instruments to Refocus; an addition of 
current results in a longer structure, and a subtraction of current from a stretched tether 
system shortens the structure. In a high k * system, the system could not be extended without 
enormous force beyond its unstretched length, and retraction would basically require 
formation flight control of the coils as they would not be constrained in their motion. 
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Stretched tethers present a risk to the spacecraft, though: if there is a power failure, even 
momentarily, the electromagnetic force keeping the tethers stretched would disappear, and 
the spring force would pull the coils back together with no electromagnetic force to slow them 
as they approach collision. This could result in damage to the spacecraft or deployed structure 
unless redundancy is built into the power subsystem for each coil to provide temporary power 
or a slow fade of current to reduce the snapback from the stretched tether. Therefore, the 
most significant trade is risk of damage versus almost instantaneous Reconfigurability. 



Force coefficient y 


Figure 32: Equilibrium length L eq * as a function of y plotted for different values of k* from 10 5 fc 0 * to 10fc 0 ’. 

For low ys, the equilibrium length does not change for all but the smallest k*. This is because 
the electromagnetic force is dwarfed by the restoring spring force in the case of larger k*s, so 
only the smallest spring forces (blue, green) are of comparable magnitude to the smaller 
electromagnetic forces and allow some stretching. As y increases, the electromagnetic force 
begins to dominate the spring force and stretch the tethers in steady state, starting with the 
smaller k*s and finally stretching the system with the strongest spring force (black) around ten 
times the baseline k*. We see a growing L eq * and difference between the different k* values 
for large ys because high electromagnetic forces dominate spring forces (which increase 
linearly as L eq * increases) and thus stretch the tethers more the smaller the k*. 
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4.3. 2. 2. 2 Varying unstretched tether length L t * 

Figure 33 shows how the equilibrium length L eq * varies as a function of y for variable tether 
length L t *, with b*, k*, and d 0 * fixed at their baseline values. The baseline value of L t * is 
pictured in teal. As per the discussion in Section 4. 3. 2. 2.1 above, we observe that the 
equilibrium length increases with higher force coefficients y and is more sensitive to the 
unstretched length of the tether at lower ys, since the spring force dominates the 
electromagnetic force there. High ys dominate the spring force, and as the tethers stretch, 

L eq * ~ L t * /L t * , or the amount the tethers stretch (A L) as a fraction of the unstretched length, 
gets smaller. For the same large y, when L eq * is large as a consequence of the large force, A L 
does not vary much with L t * , which is why the lines converge for high y. 



10 2 

Force coefficient y 


Figure 33: Equilibrium length L eq * as a function of y plotted for different values of L t * from 6 to 18 coil radii. 


Equilibrium length L as a function of force coefficient y 
for various L*and baseline b*, k*, d . * 

t ’ start 
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Varying damping coefficient b or initial distance d start * has zero impact on the steady state 
result of the system, as can be seen in Figure 34; all values of damping and initial distance 
settle out to the same steady state length L eq * in the end, regardless of the time taken to do it. 
The blue line in Figure 34 is the same as the baseline teal line in Figure 33. 


Equilibrium length L * as a function of force coefficient y 
for any b * d stan * and baseline k* 
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4.3.3 Hinged constraints in Unfold configuration 

The hinged case presented in this Phase I report examines the behavior of two coils affixed to 
panels joined with a hinge and a torsional spring and damper along a mutual edge, or the 
configuration pictured in Figure 13. The following trades are intended to provide insight into the 
selection of design variables like panel size, effective current, spring constant, rotational damping 
coefficient, and natural spring angle. 

4.3. 3.1 Moment f, versus deployment time T 

The deployment time of two hinged coils depends on the electromagnetic moment between 
them, and the damping coefficient /?* and spring constant k* of the hinge. The electromagnetic 
force depends on two parameters in the nondimensionalized equations of motion: f and the 
angle \9 X — 9 2 \ between the coils. Deployment time in the hinged case is defined as the time 
the system takes to reach the steady state. For numerical calculations, we define the steady 
state to be the case in which both the relative angular velocity and the angular acceleration of 
the two coils are arbitrarily close to zero: 


die i - 0 2 ) 


d 2 {9 i 0 2 ) 

dt 

— €i, 

dt 2 


Again, we choose: 


e v e 2 — 2 x 10 6 


The following plots investigate how the deployment time T for the Unfold configuration is a 
function of S' P*. K*' dinit ' and 9 0 . These variables encompass most of those variables from 
Table 6 that affect the motion and steady state angular separation of the hinged system 
without our plots becoming redundant with variables that are closely related to one another. 
In each of the following cases, one variable is being varied, and the others are fixed at the 
baseline values given in Table 8. 


4.3. 3. 1.1 Varying damping /? 

Figure 35 shows how the deployment time T varies as a function of f for variable damping 
coefficient /?*, with spring constant k* fixed at its baseline value. The baseline value of /?* is 
drawn in red. The longest deployment times T are for both overdamped systems (blue, black, 
and green lines) and underdamped systems. The shortest Ts are for those systems for which /?* 
for a given k * results in close to critical damping behavior (magenta and teal lines). For the 
overdamped cases (black and gold), higher fs result in the system being forced to asymptote 
more quickly towards equilibrium. In the underdamped case (blue), increasing <f simply 
increases oscillations and results in the system settling to steady state after increasingly longer 
times. The blue line experiences ranges over <f where T does not increase with f, followed by 
ranges where it does. This likely reflects a combination of several factors: an overly restrictive 
criterion for steady state deployment time (e 1 , e 2 ), the fact that f is plotted logarithmically (as it 
is in all of these figures), and the relative insensitivity of the system to the moment coefficient 
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at low values, all of which could combine and result in the presence and emphasis of numerical 
noise, especially for small f. 

This plot indicates to system designers that for a given k*, it is possible to select an ideal, close 
to critically-damping /?* (teal and magenta) value that remains ideal regardless of ^s and its 
internal variables. 

Deployment time T as a function of moment coefficient £ for various p* and baseline tc* 



Figure 35: Deployment time T as a function of f; plotted for different values of 6* from 10 5 /? 0 * to 10 fi 0 * . 
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4.3.3. 1.2 Varying spring constant k 

Figure 36 shows how the deployment time T varies as a function of f for variable spring 
constant k*, with j3* fixed at its baseline value. The baseline value of re* is not shown in this 
plot; however, the closest value is lO 009 x k 0 *, drawn in teal. The lowest k * case plotted here 
(jagged blue) takes longest to reach steady state for low moment coefficients, because the 
system is relatively underdamped with respect to the baseline /?*, and, therefore, oscillations 
are allowed to persist for a longer period of time before reaching steady state. For low k*, the 
system is more sensitive to various combinations of the baseline /?* and f, resulting in the 
jagged nature of the line. For progressively larger k* s, the deployment time is relatively 
insensitive to k* itself but increases with f and is smoother the higher k*. 

System designers will want to, for <f > .1, choose a low k* (like jagged blue) for faster 
deployment time; for f < .1 , it is advantageous to choose a k* larger than the baseline value. 

Deployment time 7~as a function of moment coefficient c, for various k* and baseline (3* 



Figure 36: Deployment time T as a function of ^ plotted for different values of k* from 10 75 k 0 * to 10 1 5 k 0 *. 


85 


Massachusetts Institute of Technology 


Space Systems Laboratory 


Electromagnetic deployment of space structures 


NASA NIAC Phase I Final Report 


4.3. 3. 1.3 Varying initial angle 9 0 and natural angle 9 init 

Figure 37 and Figure 38 show how the deployment time T varies as a function of f for variable 
starting angle 9 init with /?* and k* fixed at their baseline values. In Figure 37, the natural angle 
of the spring 9 0 — n/8, and in Figure 38, 9 0 — 3n/A. These two 9 0 values were selected to 
bracket a large range of natural angles. The baseline value of both 9 init and 9 0 is 7t/4, drawn in 
this plot in teal. Starting the coils at the same angle as the natural angle of the spring (red, or 
7r/8) results in the shortest deployment times for f < 0.1; this is because small moments do 
not stretch the spring much beyond its natural length, and therefore any disturbances from this 
9 0 are smaller and the oscillations of lower amplitude than when 9 init =£ 9 0 . For high f, 

B iri it > do results in the shortest deployment times, because the spring force is already 
working to resist the electromagnetic repulsion at the start of deployment, so that the system 
can "settle" into equilibrium relatively quickly. An exception to this observation is 9 init — 

3n/A (gold), which is so large that the spring force dominates the electromagnetic force and 
causes the system to oscillate too much. 

Deployment time Tas a function of moment coefficient q for various 0 and 9 Q = jt/8 



Figure 37: Deployment time T as a function of ^ plotted for different values of d lnlt from jt/ 32 to 3 jt/ 4 and 
0 O = n/8. 
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In Figure 38, where the natural spring angle 9 0 — 3tt/4, all of the 9 init s plotted are smaller 
than or equal to the natural spring angle, meaning that the spring force is assisting 
electromagnetic repulsion from the start of deployment. Therefore, the system will general 
exhibit oscillations that are larger than the corresponding cases in Figure 37, resulting in the 
overall deployment times observed below. We also observe that the optimal 9 init (for lowest 
deployment time) increases with increasing <f from 9 init — n/16 to 7t/2, suggesting that for 
higher moments, it is preferable to begin at an angle closer to the natural length of the spring 
but not so large that it causes the system to oscillate too much. 

From a systems standpoint and, for most applications, it is desirable for initial deployment to 
start the coils as close together as possible, to maximize the stowed-to-deployed size ratio of 
the structure. In that case, while deployment time might not be a primary concern, long 
deployment times imply large amplitude oscillations at high moments, which are not desirable. 
A current controller to reduce oscillation seems from these figures to be important for coils that 
start very close together and thus will experience high initial electromagnetic forces. One 
reason that the initial angle might be much greater than 0 is for Reconfigure or Deform 
operations, wherein the panels remain less than n apart, or not quite flat, in order to enable 
refolding or retracting of a deployed panel via electromagnetic attraction. For those 
applications, these plots are very useful. 


Deployment time Tas a function of moment coefficient for various 0 jnjt and 6 q = 3ii / 4 



Figure 38: Deployment time T as a function of f; plotted for different values of 0 lnlt from n/32 to 37 t/ 4 and 

0 O = 3n/4. 
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4.3. 3. 2 Moment % versus equilibrium angle 

As in the tethered case, the equilibrium angle 9 eq defines the final hinge geometry after 
deployment. 9 eq is independent of initial conditions, and damping and is therefore only 
dependent on the moment coefficient £ the spring constant k*, and the unstretched natural 
angle of the spring, 9 0 . 

4.3. 3. 2.1 Varying spring constant k* and natural angle of spring 9 0 

In Figure 39, we see how the equilibrium angle depends on the moment coefficient for various 
values of k* and for two different values of 9 0 , 0 and ti/4. For both values of 9 0 , the largest 
9 eq s are achieved for the smallest spring constants, indicating the inability of the springs to 
resist larger electromagnetic moments. Subsequently, the system will tend to settle to an 
equilibrium angle of 180° (blue, green, and red), where the electromagnetic moments also 
vanish. Similarly to the tethered case, smaller equilibrium angles result from larger spring 
constants. The effect of increasing the natural spring angle to ti/4 is to shift up the equilibrium 
angles from a minimum of 0 for 9 0 — 0 to a minimum of ti/4. 

Since equilibrium positions do not depend on inertia, any combination of k* and f below can be 
chosen to effect a desired equilibrium angle. System designers can select a k* and f for their 
design using this data. 


Equilibrium angle 0^ as a function of moment coefficient § for various k* and 0 Q = 0 and n / 4 




Figure 39: Equilibrium angle 0 eq as a function of moment coefficient f for various k* and 6 0 = 0 and n/4. 
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4.4 Conclusion 

In this chapter, we have made use of the models originally presented in Chapter 3 in order to 
perform high-level analysis of the interconnectivity of the system variables. This analysis was 
enabled by introducing a number of dimensionless parameters which combine many of the 
design variables and facilitate systems-level analysis by reducing the dimensionality of the 
variable space. We focused on two primary metrics: deployment time T, and the equilibrium or 
characteristic length L eq (or angle 6 eq ) of the system after deployment. The effect of varying 
design variables and their dimensionless counterparts was presented in such a way as to 
facilitate the decision-making process of a system designer. 
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5 Chapter 5: Viability 

Revisiting the initial questions of the study, we see that Chapter 4 addressed Question 1, the 
technical feasibility question: 

Question 1: Can we use electromagnetic forces generated by and acting between high- 
temperature superconductor current-carrying coils to move , unfold, and support parts of 
a spacecraft from its stowed position? 

by exploring the underlying physical concepts of electromagnetic structures and the effect of 
system parameters on deployment dynamics and steady-state positions. From Chapter 4, we 
conclude that the baseline physical principles of HTS structures are sound and that limitations 
on the use of the technology will come from mission requirements and available subsystem 
budgets. Questions 2 and 3 both involve the mission and economic viability of HTS structures in 
the space acquisition environment: 

Question 2: For which operations does this technology represent an improvement over 
existing or in-development options? 

Question 3: What new mission capabilities does this technology enable? 

and will thus be included in this chapter on the viability of HTS structures as a structural 
technology for space systems. A summary of the general advantages and disadvantages of HTS 
structures as compared to "traditional" structural options (panels, beams, motors, pyrotechnic 
bolts) at subsystem levels is included below in 
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Table 9. Red text indicates a disadvantage, or area where HTS structures are competitively 
weaker than traditional structures and mechanisms. Green text indicates an ability or property 
that HTS structures have that traditional structures do not. Black indicates a parameter that 
depends on the properties of the electromagnetic structure design. 
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Table 9: Subsystem advantages and disadvantages of HTS structures over traditional options 


Subsystem or process 

Example impacts of HTS structure choice 

Avionics/ Comms 

Potential electromagnetic interference 

Thermal 

Additional thermal control required : heatpipes + cryocoolers, 

little or no conduction 

ADCS 

Current regulation needed. 

Earth's magnetic field, 

additional momentum trading possible 

Optical Path (Payload) 

Risks to position accuracy and disturbance control 

Structures 

Reduction of mass (extra thermal mass for each coil cooling system, 
but for large structures could be an overall reduction of mass), 

vibration isolation, 

increased compaction/size ratio, 

enables reconfiguration with no additional components 

Propulsion (if 
formation flying) 

No propellant required, 
eliminates thruster plumes 

Power 

Additional power draw for both coil currents and cryocoolers 

Testing 

Difficult in lg without a complex testing rig 


There are three kinds of trades that lead to the selection of an HTS structure design. The first is 
a mission performance trade, where mission destinations, orbital parameters, and performance 
needs like reconfigurability are traded based on the high level requirements of the project. 

Since spacecraft environment is important to the use of HTSs thermally, and HTSs can enable 
some more difficult performance needs, these trades inform the next kind of trade, which 
occurs at the subsystem level. At this level, an HTS structural system is compared to other 
structural options, both traditional and "alternative" (defined below), and the impact of the 
structural choice on subsystem budgets is evaluated. At the final, component level, trades are 
performed amongst various HTS designs, such as those in Chapter 4. 
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In this chapter, we will discuss in a little more detail the comparison between HTS structures 
and other "alternative" structural technologies like inflatables and tensegrity structures. We 
will also explore some of the design trades that affect the "internal" variables, or the 
constituent variables of the dimensionless parameters plotted in Chapter 4, as well as the other 
variables included in Table 6. 

5.1 Alternative structural technology trades 

Currently, structural designers reduce structural mass by using lightweight materials, such as 
aluminum frames, and bydesigning for a high strength-to-weight ratio using isogrid panels and 
structural optimization. In order for HTS structures to find a foothold amidst these traditional 
structural options and other innovative, alternative solutions focusing on mass reduction, they 
should present a strong advantage over other techniques in important metrics like mass, 
compaction ratio, and thermal and vibrational isolation. Then, despite the additional thermal 
and power needs of HTS structures, they could still be a preferable selection. 

In 
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Table 10, we summarize important aspects of several leading alternative structural 
technologies. "Alternative" is taken to mean a departure from purely rigid, metallic structural 
elements in this discussion, such as composite materials (both rigid and shape memory), 
tensegrity structures, and inflatable membranes. 


94 


Massachusetts Institute of Technology 


Space Systems Laboratory 


Electromagnetic deployment of space structures 


NASA NIAC Phase I Final Report 


Table 10: Description of alternative structural technologies 


Alternative 

Description 

Contributors/ 

Disadvantages w.r.t HTS 

Structural 

Technology 


Projects 

deployed structures 



Composites 


Lightweight and strong • SpaceShipOne 

Customizably anisotropic 
Maintains shape • ELV payload 

Good for enclosures fairings 

Can reenter atmosphere • Pellegrino 

Low coeff. of thermal (Caltech) 

expansion 


• No thermal or vibrational 
isolation 

• No staged deployment 

• Mostly traditional 
mechanisms for 
deployment 

• No reconfigurability 


Tensegrity • Very lightweight and stiff 

• Keeps beams under 
compression with 
tensioned cables 

• "Snap out" non- 
reversible deployment 



Knight, B. 

Fuller, B. 
Skelton, R. 
Snelson, K. 


• No thermal or vibrational 
isolation 

• Polygonal shapes - may not 
be advantageous for 
mirrors or lenses 

• No reconfigurability 

• Size limited by beam and 
tether lengths 

• No staged deployment 



• Very lightweight material 

• Easy to implement round 
shapes 

• Inflates a membranous 
shell with gas 

• Compact stowage 


L 


Echo 1 
Explorer 9 
L'garde ( 

Antenna 

nt, 

decoys) 
Adherent 
Cadogan (ILC 
Dover) 


• Leakage overtime (UV 
curing might help) 

• One rip renders useless 

• No reconfigurability 

• No staged deployment 

• Compressive loads require 
pressure or curing 

• Must carry pressurized gas 
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5.1.1 Comparison of HTS structures to alternative structural technologies 

Table 11 qualitatively lists key trades between HTS structures and other structural support and 
actuation technologies. The Other Technology column lists the technology against which HTS 
structures are being compared. The second and third columns show the desirable attributes of 
the Other Technology versus HTS Structures respectively. 


Table 11: Comparison of HTS structures to alternative structural technologies 


Other Technology 

Other Tech Attribute 
vs. 

HTS Structures Attribute 

Inflatables 

Low power and 
potentially lower mass 

Longer lifetime 

Tensegrity 

Low power and no 
thermal control 

Controllability of deployment and potential 
reversibility 

Composites 

Low power and no 
thermal control 

Shape change and compaction/size ratio 

Pyrotechnics 

Reliable impulse 

Deployment control and potential reversibility 

Piezoelectrics 

Precision 

Range of motion 

Spring-loaded 
booms and hinges 

Lower power and 
potentially lower mass 

Controllability of deployment and potential 
reversal 

Traditional 

motorized 

actuators 

Heritage 

Low mass and thermal/vibrational isolation 


Table 12 depicts a qualitative comparison of electromagnetic support and the three structural 
support technologies listed above (inflatables, tensegrity structures, and composites) across a 
variety of characteristics of interest for spacecraft designers. Values of these characteristics (in 
this table, yes, no, or other) that are generally considered good for the specific metric by 
spacecraft designers are green, values that are dependent on the specifics of the 
implementation of the technology and could be either good or not are yellow, and those values 
that are considered less desirable are red. All characteristics are treated equally in this table, 
but in some applications, some attributes may be valued more highly than others. The 
takeaway from Table 12 is that HTS structures compare favorably to other structural techniques 
in most performance metrics associated with spacecraft structures, but that the HTS power and 
thermal requirements are significantly less desirable. Further development of HTS structures 
should focus on mitigating the risk that the power and thermal metrics present to successful 
technology integration. 
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Table 12: Qualitative comparison of deployable structure options 


not 


CHARACTERISTICS OF 
INTEREST 



Lightweight? 

Yes 

Yes 

Yes 

Yes 

Strong? 

Yes 

Yes 

If cured/rigid ized 

Yes 

High compaction ratio? 

Booms under 
development 

Yes 

Yes 

Yes 

Beams limited by payload fairing 

Only if solid panels 

Yes 

No 

No 

Traditional deployment mechanisms 
required (mass) 

Often 

Yes, high 
loads/shock 
otherwise 

Yes 

No 

Durable? 

Yes 

Yes 

No 

Yes 

Vibrational isolation? 

No 

No 

No 

Yes, in some configs 

Thermal isolation? 

No 

No 

No 

Yes 

Low coeff. of thermal expansion? 

Yes 

If made with 
composite bars 

Only if shape memory 
polymer (SMP) 
composites 

Yes 

Thermal control required? 

No 

No 

No 

Yes 

Power required for deployment? 

No 

Depends on 
method 

No 

Yes 

Staged deployment/ assembly? 

No 

No 

No 

Yes 

Can be non-polygonal? 

Yes 

No 

Yes 

Yes 

Reconfigurable w/out mechanisms? 

No 

No 

No 

Yes 


5.2 Internal variable trades 

The variables internal to the nondimensionalized coefficients used in Chapters 3 and 4 
encompass other variables in the design vector for an HTS structure, including mass (of coils m, 
of attached components like m pane i, or of the total system m tot ), power draw (of each coil, 
taken as P co u, or of the thermal subsystem, taken as P t hermai)> number of turns in the coil n, 
diameter of coil D, stowed structure length L s , the damping coefficient b or /?, the spring 
constant k or k, the temperature that a coil is maintained at T e , the magnetic field B, and the 
different forces and torques. 
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Table 13 displays several of these system trades: the first column names an increase or 
decrease in one variable, and the second column reflects the reactions of other system 
variables to that change. Red text is a change that is generally not in the direction expressed in 
an objective vector for spacecraft (an example general one shown below), such as an increase 
in mass or power. Green text is a change that is in the direction specified in the objective 
vector, such as a minimization of mass or power. Black text labels variables that can optionally 
be changed along with the first column change and which might be good or bad or neither. In 
the case of increasing dimensionless spring or damping constants, for example, selection of 
either of these variables is likely not due to maximizing or minimizing these variables 
individually but rather through the maximization or minimization of another variable related to 
both constants from the plotted relationships in Chapter 4, such as equilibrium length. 


r 




EXAMPLE 

OBJECTIVE 

VECTOR 



min(total mass) 
min(power req'd) 
min(deployment time) 
min(thermal mass) 
min(mass/length or area) 
min(deployment steps) 
max(size ratio) 
min(current switching) 
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Table 13: System variable trades 


Trade 


Temperature 

of coil T e si/ 

in order to 

increase 

critical 

current 

density 


Effects on other variables 


/ f" Thermal power draw 'T' 

I s Magnetic field to mass ratio I s 
I s Coil power draw / f s 


Variables 

involved 


Pcoih 
P thermal' 

T e , B, m 


Diagram 


S Cold ^ 


Number of 
turns n in coil 
^ to increase 
effective 
current 


'f- Magnetic field from higher 
effective current 'T 

I s Mass of coil I s 

si. Current through each wire si/ 

Power/mass ratio 


B , m coi i, 

Pcoib P/ 


eff 


si/ Mass of coil si/ 

I s Magnetic field '!' 

1" Force 'V 

Dimensionless damping 
Dimensionless spring constant 




i i 

Less turns in coil 




More turns in coil 


Radius of coil 
R si/ to 
increase 
magnetic 
force 


R, m co u, 
B, F, k or 
k, b or /? 



Larger coil 


Smaller coil 


# of coils N ^ 
to reduce 
current in 
coils or 
increase 
characteristic 
length 


/ f' Mass of system I s 
'l' Thermal power draw 'T' 
and 

- Linear structural density - 

'V Characteristic length I s 
‘f- Coil power draw I s 
or 

'I s Linear structural density 'T' 
- Characteristic length- 

si/ Coil power draw si/ 


EP-tot’ L e qi 
Pcoih 
P thermal 


More coils 
Same current 
OR 

Longer L eq 


Less coils 
Smaller current 
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Amount of I s Coil power draw / f s 
current / ^ 'T' Magnetic field 'T' 

I' Force 'V 

't' Characteristic length / \' 

4/ Size ratio ^ 

Dimensionless damping 
Dimensionless spring constant 


Pcoih L e qi 
B, F, k or 
k, b or /? 


Same L eq 
Smaller current 




Some more detailed discussion of the trades in Table 13 is included below: 

• Temperature of coil T e ^ 

o This trade makes reference to the relationship between the critical current density 
that an HTS wire can support and its temperature. As temperature decreases below 
T c , the critical current density in the wire also increases, so at colder temperatures, 
more current can be run through the same coil, and thus a higher magnetic field and 
force can be generated at the cost of more thermal power (to cool the coil to a lower 
temperature) and more coil power (to generate the higher current) 

• Number of turns in coil n I s 

o Increasing the number of turns in the coil increases the effective current in the coil to 

hff = 

and therefore increases y, increasing the magnetic force. However, when current is 
flowing through multiple turns of a coil, the magnetic field of one turn exerts a 
magnetic force on the current in another turn, forcing it towards one side of the wire 
and effectively reducing the size channel through which the current can flow. The 
wire's critical current density limits how much current can flow through any portion 
of the wire's cross section, so once that limit is reached, the total current cannot 
increase; therefore, additional turns of wire reduce / and result in a loss term. 

The same change in power per unit mass value can result from an increase in mass, 
with power being held constant, or from a decrease in power, with the mass being 
held constant. The latter would indicate a decrease in effective current unless n is 
being increased and R is being decreased to result in the same mass but a higher 
possible effective current, as observable in the relationship: 

Peon IV 
m coil n2npR 

where V is the voltage across the coil and is assumed constant. 
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• Radius of coil /? sj/ 

o We recall that the mass of the coil may be written: 


m — 2nRp 

p is the linear density of the coil, given by: 


p = p v nA c 


and p v is the volume density of the coil material, n the number of turns in the coil, 
and A c the cross sectional area of the coil. We may then rewrite the dimensionless 
force parameter yas follows: 

^ 8n 2 pR 2 


Hence, we see that: 


Y 


1 



so that decreasing the radius of the coil increases the magnetic force. 


While it may seem like making R arbitrarily small is beneficial due to the increase in 
the magnetic field through y, we must keep in mind that changing R also has an effect 
on the dimensionless damping and spring coefficients, b* and k*, respectively: 


2npR 

kT 2 
u* — 

2npR 

Thus, if the designer has a limited choice of actual materials, and hence b and k, in 
order to implement a desired dynamical and structural configuration, R cannot be 
made arbitrarily large or small and must in fact be chosen so that the interplay of 
damping and elastic forces give the desired deployment dynamics and equilibrium 
configuration as per the analysis in Chapter 4. 


• # of coils N't' 

o The number of coils in the system N can be increased for one of two primary 

purposes: the increasing of the characteristic length of the structure with the same 
current per coil, or the reduction of current through each coil with the same 
characteristic length. The first purpose results in no increase in the linear structural 
density, or mass per unit length, whereas the second purpose increases the linear 
structural density as there are now more coils for the same characteristic length. This 
dual choice is reflected in the bracketed terms in Table 13. Either choice results in an 
increase in the mass of the system and in the power required for cooling of the coils. 
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• Amount of current / I s 

o For a current density below the critical value for the wire, increases in current result 
in an increase in the magnetic field and thus electromagnetic force. For a given 
damping and spring constant, the characteristic length or angle will then increase as a 
result of the increased electromagnetic force. Thus, the size ratio (stowed-to- 
deployed) of the system will decrease but the coil power draw will increase due to 
the increased current. 

Now that we have set up the basic constraints on coil motion, in future work we can now 
expand on these constraints to include multi-operation deployments, complex tether 
configurations and multi-hinge cases that enable Reconfiguration or temporary Deformation of 
a paneled structure. These studies, to occur in the course of our Phase II work, will examine the 
feasibility of "new mission capabilities" as motivated by Question 3. 
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5.3 Phase II study plans 

Our proposed Phase II effort will, in the two additional years of study time beyond the end of 
this Phase I study, expand upon the primary questions of our Phase I study, develop HTS 
structures further towards the attainment of TRL 3, and reduce the risks associated with new 
technology development. Our Phase II effort includes Prof. Raymond Sedwick, Director of the 
University of Maryland's Space Power and Propulsion Lab, as a co-investigator focused on the 
development and testing of a flexible cryogenic heatpipe as a key enabling technology for a 
number of HTS structural applications that require a flexible HTS coil (such as those in the 
Expand configuration). We also add collaborators from the NASA Goddard Space Flight Center 
(GSFC) Mission Engineering and Systems Analysis division, including Dr. Carl Adams, assistant 
chief of technology, to the team to support technology trade studies and integration feasibility 
assessment and analysis for HTS structures in a large, complex, and fairing-limited precision 
structure such as the James Webb Space Telescope (JWST). 

In our Phase I work, we have identified a number of potential barriers to the feasibility of HTSs 
as a replacement for traditional structures. We have sorted them into three categories of 
feasibility barriers to be addressed in Phase II in order of decreasing impact: Power and 
Thermal, Dynamics and Control, and Implementation. The four impact categories with which 
each concern is associated: Cost (dollars, mass and power). Spacecraft Performance, 
Development Time, and Key Enabling Technologies. The three key enabling technologies that 
our technology development effort focuses on are rigid deploying HTS coils, flexible deploying 
HTS coils, and flexible cryogenic heatpipes. All of these key technologies require experimental, 
modeling, and feasibility work in our Phase II effort. 

The power and thermal subsystems are those most adversely affected by the use of HTSs, as 
described earlier in this report. Because a large portion of the power required by HTS 
structures is consumed by thermal equipment, these two subsystems are inextricably linked in 
our study of feasibility issues. The power and thermal barriers we will discuss in this section 
are: 

• Cost: Complex power and thermal control, generation and management, 

• Cost: Thermal vacuum chamber testing expenses 

• Performance: Power for HTS operations needing to be taken from other subsystems 

• Key technology: Low-power rigid and flexible cryogenic heatpipes 

• Performance: Lower power cryocoolers 

Coil dynamics and control with coil boundary constraints, though, introduce performance 
feasibility issues that must be addressed in order for any technology advancement to take place 
and include: 

• Performance: Coil dynamics 

• Performance: Stability of configuration 

• Performance: Position accuracy and knowledge 

• Performance: Disturbances (Cryocooler vibration, solar pressure) 
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Implementation feasibility concerns center around the competitiveness of HTS structures with 
other technologies, either in cost, key performance metrics, or the ability of less complex HTS 
structures to be competitive enough of a structural option to find purchase on missions in order 
to gain flight heritage to then motivate the research and development of more complex HTS 
structures. 

• Cost: Cost may be uncompetitive with other structural options 

• Performance: Performance may be uncompetitive in multi-subsystem tradespace 

• Performance: Potential EMI, negatively affecting other subsystems on board or vice 
versa 

• Development Time: Technology roadmap may be too linear to allow missions to use less 
complex HTS structures as flight hardware prior to completion 

Two more years to study the implementation of HTS structures under a NIAC Phase II allows 
us to: 

• Reduce risk in areas of uncertainty with regard to feasibility issues above 

• Assess the mission-general and mission-specific (JWST) impact of many of the feasibility 
issues discussed in this section with simulation and experimentation 

• Collaborate with NASA to fully understand the needs of the customer in order to lay out 
a reasonable and accomplishable roadmap towards integration 

• Advance key technologies necessary for full utilization of HTS structural configurations 

• Explore the impact of HTS structures on other spacecraft subsystems in greater detail, 
elaborating on the trades suggested in Phase I with quantitative comparisons 

5.4 Summary and conclusions 

In this Phase I report, we have conducted an initial feasibility study for the use of high- 
temperature superconducting coils as deployment and support elements in spacecraft 
structures. We have also introduced system performance trades that will be further quantified 
in the Phase II work. 

In Chapter 1, we introduced the notion of high-temperature superconductors and the ways 
they could potentially be useful in spacecraft structures. We discussed the five primary 
benefits that HTS structures can offer to the aerospace community: 

1. Reduced mass 

2. Larger structures with same launch vehicles 

3. Vibration- and thermally-isolated structures 

4. Staged deployment, in-space assembly with fewer operations, and partial system 
replacements 

5. Reconfiguration of structures after deployment 

and introduced the three primary research questions relevant to this study. 
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In Chapter 2, we motivated the work of Chapter 3 by discussing the physical principles 
underlying superconductivity and how electromagnetic interactions between conducting coils 
can result in their motion. We also discussed previous work in the use of electromagnetic 
forces for relative motion, positioning and deployment as well as the state of the industry for 
superconducting wire construction. 

In Chapter 3, we developed the equations of motion for coils responding to electromagnetic 
forces while under the influence of constraining elements (i.e. tethers and hinged panels) and 
validated our numerical models against known analytical solutions. 

In Chapter 4, we modified our equations of motion for trade analysis by nondimensionalizing 
them and introducing several dimensionless parameters that encode many of our design 
variables. We presented the results of numerous simulations using the numerical models 
introduced in Chapter 3, detailing the response of tethered and hinged coil systems for various 
dimensionless parameters and, hence, combinations of the original design variables. 

In Chapter 5, we discussed on the basis of the results from Chapter 4 the viability of HTS 
structures in the context of trade analyses. Trades were described at the mission level, the 
structural subsystem level, and the component level. 

On the basis of the work performed in this Phase I study, we have come to the following 
conclusions: 

1) The fundamental physics support the basic concept of using electromagnets for 
deployment. 

2) With existing technology and capabilities, the use of HTS coils for the deployment 
activities investigated herein is achievable, though with significant work to be done in 
dynamic controllability, system optimization, sizing and scaling, thermal control, and 
mission integration analysis. Our Phase II work will address many of these topics as well 
as economic viability. 

By the end of Phase II, we aim to have completed and tested an experimental flexible cryogenic 
heatpipe, a technology integration case study of JWST, laboratory experimental demonstrations 
of HTS coils unfolding panels and separating while superconducting, and models calibrated both 
by experiments and analytical solutions. These are the first steps along the HTS technology 
demonstration path, and those steps combined with the foundation laid in our Phase I report 
and the feasibility and implementation analysis to be conducted in Phase II will bolster the 
credibility and prepare the technology for further TRL advancement. 
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